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PBEFACE 


The present report contains results of an investigation of critical 
problems related to satellite systems for long-term earth energy budget 
(E£B) observations, performed under Langley Research Center Contract No. 
NASI-11871 for the National Aeronautics and Space Administration. 

This research study, performed by Drexel University, is one part of 
a much larger effort by several institutions, including Colorado State 
University, The University of Wisconsin, Virginia Polytechnic Institute 
and State University, Link Temco Vought as well as cognizant personnel at 
NASA Langley Research Center. This team is studying satellite systems 
for performing long-term EEB measurements over geographical regions, hemi- 
spheres, and the entire earth for periods of 10 to 30 years. A major 
portion of the total effort is responsive to the AAFE proposal, and the 
proposed LZEEBE system (reference 1) which employs three balloon 
radiometers . 

The decision was made to expand the scope of the total effort beyond 
that envisioned in the AFFE proposal. This boradened scope Includes Phase 
A type of efforts concerning other concepts of performing EEB observations 
besides the balloon system. For example, systems employing spinning plate 
radiometers and/or scanning radiometers could be developed for long-term 
space application. Regardless of the geometric characteristic of the 
observational system, the problems of data analysis and interpretation are 
similar for all wide-angle systems with only an adjustment required for 
viewing geometry. 

The current investigation was performed during the twelve month period 
1 January 1974 through 31 December 1974. This period is denoted as Phase III 
in subject contract. The express purpose of the investigation is outlined 
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in the following objectives, contained in the statement of work. 


1. Accuracy assessment of Sampling Bias by Candidate 
Satellite Systems for EBB Observations 

2. Development of Procedures for Analysis and Inter- 
pretation of EEB Observations by Wide-field, Broad- 
band Detectors 

3. Accuracy Assessment of Procedures for Analysis and 
Interpretation of EEB Observations 

4. Specification of Spectral Absorbing Characteristics 
of Broad-band Detectors and Calibration Requirements 
for Accurate Observation of the EEB 

5. Variance Spectrum Analysis of EEB Observations from 
ESSA VII Satellite 

6. Professional Support and Advisory Effort as Co-Principal 
Investigator of LZEEBE 

This report is the final report for the phase III effort under NASA 
Contract No. NASl-11871. Other reports that are related to the overall 
effort are: "An Investigation of ESSA VII Radiation Data for Use in Long- 

term Earth Energy Experiments," published as NASA CR-132623; "Our Con- 
taminated Atmosphere - The Danger of Climate Change," published as NASA 
CR-132625 and "Steady-state Solution to the Conduction Problem of a Spherical 
Balloon Radiometers," published as NASA CR-132624. 

Gratitude is extended to several NASA/LaRC personnel for their encou- 
ragement, Interest, stimulating discussions and suggestions provided during 
the present investigation. Among these scientific personnel are Included: 
Messrs. George Sweet (technical monitor), Charles Woemer, Jack Cooper, 

Dr. Louis Smith and other members of the LaRC team. 

Frederick B. House, Project Director 
Associate Professor of Physics and 
Atmospheric Science 
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INTRODUCTION 


The study of the total energy budget of the planet earth and its atmo- 
sphere is based on analysis of the exchange of radiant energy between the earth- 
atmosphere (E-A) system and space. Radiometers on earth-orbiting satellites 
can effectively measure this exchange of radiant energy at observation points 
in space external to the E-A system. The current effort is concerned with the 
analysis and Interpretation of observations by wldeangle, spherical and flat 
radiometers with regard to the general problem of the earth energy Budget (EEB); 
but in particular, it is concerned with the problem of determining the energy 
budget of regions smaller than the field of view (FOV) of these radiometers. 
Before considering specific reasons or motivations for conducting the present 
research effort, the background to the overall problem will be discussed. 

Background 


The energy budget at a specific time t for a given region of the E-A 
system's surface can be described by the following expression which shows the 
relationship of the net flux Q to the three fluxes H , W , and W . 

S 1l g 


Q = H - (W^ + W^) 
s re 


(1-1) 


The meaning of the three fluxes is shown schematically in figure 1-1. The 
symbols Q, W^, and are defined in the List of Symbols in the front matter 
of this report. Introducing the concept of albedo A = W /H , which is also 

1C s 

depicted in Figure 1-1, . one .can.; rewrite (1-rl). as 


Q = (1 - A) - 

s e 


(1-2) 
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Figure 1-1. Radiant power balance on an earth surface element. 
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A time average of the aet flux Q over any convenient time scale 
can be readily computed from- the data available for any of the E-A system re- 
.glons that has been under study. 

The net flux Q defined by equation (1-1) can represent the time average 
conditions for a broad range of spatial dimensions. These can range from the 
dimension S'-’Of a G ARP id, about 5^ x 5 ° great circle arc, to the dlmenslons.:^^ 
of the entire globe. Different observational satellite systems must be 
anployed to perform the required measurements over this broad range of spatial 
dimensions. For example, a scanning radiometer would be needed to observe • 
the components of net flux for a GARF grid. On the other hand, a wide-angle 
radiometer can best fulfill the sampling requirements for the entire globe. 

In terms of conventional meteorological dimensions a GARP grid maybe classed 
as a mesoscale feature. Wide-angle radiometers become more effective obser- 
vational tools than scanning (narrow angle) radiometers for some meteorological 
features between the mesoscale and global dimensions. It Is suggested here 
that the synoptic scale meteorological feature may be the cross over point 
between the usefulness of scanning and wide-angle radiometers. 

The current effort focuses on the problem of Interpreting wide-angle 

observations In terms of the components of the net flux for regions whose 

dimensions are synoptic scale and larger, l.e., spatial dimensions whose areas 
6 2 

are 1-5 x 10 km to the entire area of the earth. It should be emphasized that 
this problem of data Interpretation Is common to both spherical and plate 
radiometers. 

The Important point to be noted here Is that W^, and W^, can be 
determined from observational data of satellites orbiting the earth. It follows 
then, that Q can be computed for any region for which satellite data are 
available. Consequently, determination of the fluxes W^ and W^ mentioned above 
for regions smaller than the FOV of the satellite becomes Important, and 
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techniques to accomplish this task have been sought for some time. 

At Drexel University, two techniques involving matrix inversions were 

^developed or ^computing- the* two-fl\ixes- under- dis cus sion-for-^regions^^smaiier^ 

than the FOV. The motivation behind the efforts to determine W and W for. 

re 

surface areas smaller : than the FOV will now be discussed. 

Motivation 

One can think of two reasons for developing techniques to Interpret 
wide-angle observations for synoptic scale regions. These reasons are: 

(e) Knowledge of the radiant energy budget over several rareas 

(such as eastern continental U.S.A. , the polar caps, the 

Northwestern part of Africa, and the Caribbean - Gulf of 

Mexico region) possessing meteorological significance in the 

dynamical ahalysls of the atmosphere are of practical and 

scientific value. These areas are smaller than the FOV of 

wide-angle 'radiometers on satellites at orbital altitudes 
commonly used. 

(b) Several of the regions of Interest in studies of the earth's 
energy budget may have dimensions equal to, or larger than, 
the FOV of the radiometer (e.g . subtropical oceans). However, 
often only a portion of the FOV covers part of the area of 
Interest during an observation; and it is, therefore, impos- 
sible to determine from this measurement alone what fraction 
of the power measured originated at the area of concern. 

The above arguments provide sufficient reason for endeavoring to develop 
the types of techniques sought. The essentials of the two techniques developed 
at Drexel follow. 
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Essentials of Techniques 

In one technique, the number of observations m equals the number n of 
regions under study, and a unique solution is obtained for each region. The 
resulting system of n simultaneous equations Is solved by operating on the 
observed power column matrix with the inverse of the configuration factor 
matrix. If all the observations are made at one time during a single pass of 
the satellite, the solutions are termed INSTANTANEOUS. For this reason, the 
technique Is called Instantaneous/Inversion Technique. 

In the other technique, the number of observations m is larger than the 
number n of regions being observed. The n simultaneous equations required 
to seek a solution Is obtained by using an extension of the method of least 
squares to determine the surface which best fits the data. For this reason, 
this technique is termed the Best Fit/Inversion technique. Again here, the 
system of n simultaneous equations is solved by the use of a matrix inver- 
sion subroutine. 

The values of obtained by both techniques from error-free observations 
were totally acceptable. However, when the measurements Included uncertain- 
ties, not all the instantaneous values of computed by the first technique 
were acceptable. Nevertheless, a prediction scheme was developed which 
forecasts the quality of the instantaneous values to be determined by 
the first technique. In this scheme, the square matrix of the coefficients 
of the n simultaneous equations found by the first technique was used to pre- 
dict the acceptability of the instantaneous values of obtained for each of the 
regions. The computer program calculates the elements of the matrix and then 
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proceeds to analyze the matrix In order to predict the quality o£ the data 
to be obtained. The actual results computed were always compared with the 
pertinent tolerances listed In TABLE 7, Accuracy Requirements^ for. Radiation 
Budget Components r In order to determine the acceptability of the results. 

In all cases, the predictions were In agreement with the results of these 
comparisons. The results obtained by both techniques are now presented. 

Accuracy of Results 

Two- types ^of ‘errors ' were selected to perturb the power measurements made 

2 

by the radiometers. These were, systematic errors (0.3, 0.6, and 0.9 watta/m ) 

and gausslan random errors (the gausslan distribution had a sigma value a » 0.5 
2 

watts/m ). From TABLE 4-7, the pertinent "desired" and "minimum useful" 

2 2 

accuracy requirements for are £3 watts/m and £15 watts/m , respectively. 

To evaluate the Instantaneous/Inversion Technique, ten sets of six observations 

each were Included in the analysis. The predictions and the uncertalntltles 

In the values of obtained by using a combination of the gausslan errors 

2 

mentioned and 0.9 watts/m ;.systematlc errors were as follows (refer to 
TABLE 4-8) 



Wej^ 

We^ 

We 2 

We^ 

We^ 

Weg 

Predictions , sphere 
Prediction, plate 

Poor 

Poor 

Poor 

Poor 

Accept 

Accept 

Accept 

Accept 

Poor 

Reject 

Accept 

Accept 

Sphere 

15.2 

20.1 

6.9 

3.1 

17.4 

5.8 

Plate 

15.4 

21.9 

9.1 

4.2 

33.8 

9.3 

As It can be seen 

from these 

results 

, only the 

values 

of We 2 and 

We^ 


2 

exceed significantly the minimum useful accuracy requirements (±15 watts/m ) 
given for both satellites in TABLE 4-7. However, it Is noticed that consis- 
tently the plate exhibits larger uncertalntltles than the sphere. It should' • 
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be remarked here that different satellite orbits that Include larger segments 
within their FOV 's of regions one and five would serve to provide values 
with- lesser errors for these regions. 

In order to test the Best FitVinversion Technique, thrity-six obser- 
vations of six regions were Included. The values of obtained when the 
power measurements were assumed to exhibit the same combination of errors 
used above for the Instantaneous/Inversion Technique were as shown below. 

The six values used as a standard for comparing the values of retrieved 
were the averages of the given values computed for each of the six-regions. 




CM 

> 

We^ 

We, 

4 

We^ 


Average (Standard) 

179.4 

200.4 

220.4 

239.8 

259.6 

280.4 

Results 

179.3 

200.6 

222.2 

243.0 

257.4 

,;281.5 


As It can be seen from these data, the highest uncertainty In the value 

2 

of W^.ls exhibited by We^, which Is 3.2 watts/m . However, even this error 
is below the minimum useful accuracy requirement shown in TABLE 4 t7. 

BASIC IDEAS 

The basic Ideas and concepts required for deriving the expression for 
the total radiant power P (yhtts) intercepted by a radiometer orbiting the 
earth are discussed In the following subsections. 

4 

Shape or Configuration Factor 


This factor F which Is dimensionless, appears often in radiative 


transfer literature and is defined In the following expression (reference 3) 


P = FAW ‘ --- (2-IT 


where 

P > radiant power (watts) Intercepted by the radiometer. 

A > characteristic area of the radiometer. This area Is A *• 
ira2 for both a sphere of radius a and a flat circular 
plate of radius a. 

2 

W B radiant flux (watts/'m ) per unit time, or radiant power 
per imlt area, which Is either emitted (radiant emlttance 
W ) or reflected (radiant reflectance W ) by the earth- 
atmosphere (E-A) system. ^ 

F • shape, or configuration factor (dimensionless). This 
factor represents that fraction of an observed area fltix 
W that Is Intercepted by the radiometer per unit of 
characteristic area. 

The power absorbed by the radiometer Is given by, 


P =-ct P = a FAW 


( 2 - 2 ) 


where 


a » absorptance of the radiometer (dimensionless) . It is 

the ratio of the power absorbed to the power intercepted. 


Radiometer Characteristics 


Two types of radiometer will be treated In this report: (a) spherical, 

and (b) horizontal, flat, circular. 

2 

The characteristic area A •* ira of both of the above radiometers will be 
2 

assumed to be A > Im . This assumption serves to simplify the expressions 
without affecting the physical significance of the results. 

Both satellites will be assumed to be blackbodies l>.e.<, the spectral 
emlsslvlty and absorptivity are assumed to be unity 
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1 


(2-3) 




-Therefore, all the radiant power P. intercepted by either radiometer (sphere^ 
or plate) is totally absorbed, and one has 


P a p 


(2-4) 


Hence, the power absorbed by the radiometer can be written as 


P = FW 


(2-5) 


Radiance and Radiant Emlttance 

These two quantities are of great Importance in the discussion of long- 
wave radiation (LWR) emitted by the E-A system. One is Interested in the 

2 

expression that relates the radiance N(6,t^) (watts/m - sr) and the radiant 

2 2 

emlttance (watts/m ) of an area element dA(m ) whose centroid is located 
at longitude! wd latitude .The radiance is the radiant power per unit area 
emitted b y dA within an element of solid angle dfi. The radiant emlttance is 
the total radiant power per unit area (or radiant flux per unit area) emitted 
by dA into a 2ir st^adlans solid angle. The expression relating the instan- 
taneous values of and N for a given dA (A,<f>) is . 

W (!,<j>,t) =■ dij; N(9,i|^; X,(i),t)sin 9 cos. 9 d9 (2-6) 

e 0 0' 

Where the radiance is being considered as a fiinction of the zenith angle 6 
and the azimuthal angle which define the direction of N(6, X, (^, t). 

This quantity and ‘are functions of the position of the area observed as 
well as of the specific time at which the observation is made. In order to 
simplify the notation, the dependence of N on X, and t will not be shown. 


2-3 



A reasonable approach in Che Interpretation of LWR Is to assume Chat Che 
radiance is a function of Che zenith angle 6 only (reference ^). Then, one 
can write, 


N(9,i|;) = N(9) = N^f(9) 


(2^7) 


where » Zenith radiance, l.e., the radiance in the zenith 

direction, 9 =» Oo. 

f(9) ■ limb darkening function (LDF) which gives a measure 
of the anisotropy of the radiation field-. 


Hence, by substituting equation (2-7) into (2-6) one obtains 


W = 2tt f(9) sin 9 cos 9 d9 


(2-8) 


or 


= 


W 


2ir f(0) sin 9 cos 9 d9 

0 


(2-9) 


The integral appearing in (2-8) and (2-9) is used often and hence, it is 
concenienc to define it as. 


I (f) = f(9) 0 cos 9 d9 

0 


(2-10) 


Then, C2-8) and (2-9) can be rewritten as 
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W = 2 IT I (f) 
e ^ ' 

and 


For an ISOTROPIC radiation field, f(6) = 1, 1(f) ® Ht a^d the Isotropic 

iso 

radiance N Is then 

jjiso^_e (2-12) 

Figure 2-1 shows a satellite S at a height H above a spherical E-A system 
of radius R. The geometric nomenclature to be used in this report Is presented 
In this figure for future reference. A flat earth Is depicted tangential to 
the spherical system at the sub-satellite point (SSP). The concept of a flat 
earth will be used to make certain calculation simplifications and It will 
be thoroughly discussed In the subsection entitled "Flat Earth Application". 
From this figure one can readily write. 

6 = Y + a 

R sin 0 = (R + H) sin a (2-13) 




2ir I (f) 


W 


( 2 - 11 ) 


R sin Y = r sin a 


(2-14) 


r^ = R^ + (R + H)^-2R(R+H)cosY 


(2-15) 
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Irradlance, Radlancf v.Kadiant Reflectance, and Albedo 

These four quantities are of great importance when the solar shortwave 

radiation (SWR) reflected by the earth is treated. One of the expressions of 

2 

interest is the relation between the radiance (watts/ (® - sr)) and the 

2 2 

radiant reflectance (watts/m ) of an area element dA(m ) whose centroid 
is located at longitude X and latitude (^. The radiance in this case is the 
radiant power per unit area reflected by dA within an element of solid 
angle dO. The radiant reflectance is the total radiant power per unit area 
(or radiant flux per unit area) reflected by dA into a 2ir steradians solid 
angle . 

Before writing the expressions connecting and it will prove 

2 

helpful to define the solar Irradlance H (watts/m ) . This is the amount of 
radiant power per unit area (or radiant flux per unit area) impinging upon 
dA(X,(j)) from all directions contained within a 2ir steradians solid angle. 
Hence, the definitions of W^, H, and are very similar, except that they 
refer to emitted, incident, and reflected radiations, respectlvley . Then, 
one can write, using the nomenclature of reference 3. 

H(L, 5, X, (J), t)=SL cos C wN 


where 


2 

S » Solar constant (S » 1353±21 watts/m » 1.940±.03 cal/ 
(cm2-min) , see reference 5) . 

—2 2 — 

L " d /d .d and d are the mean and true sun-to-earth distances 
respectively. 

a Zenith angle of the sun at dA (X,4>). 

X B Longitude of the observed area element dA. 

4i B Latitude of the observed area element dA. 
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t ■ specific time at which observation was made. 


The instantaneous albedo ACC»^>4>>t) for the area element dA (X, (^) is . 
defined by - 


A(S, X, (fr, t)=W^(L, C, X, t)/H(L,(;,X,(|),t) 


(2-17) 


In order to simplify this expression, one can assume that L does not change 
significantly during the time interval of the observations. Then, for 
a given L, the dependence of and H on L does not have to be shown. In 
a similar way, X, and t can.be dropped from the. above expressions, and 


A(C) = (?)/H(?) 

From the definitions given for and at the beginning of this 
subsection, one can write the following expression relating these two 
quantities, which is similar to equation (2-8). 


W^(C)= N (9,ij;,?)sin 0 cos 0d6 

to 0 t 

Dividing this expression through by H(^), one has 


(2-19) 


A(?)= r(5) = dtp p (0 ?)sin 0 cos 0d0 (2-20) 

0 o 

Where , 

W . . 

A(0 = r(?) = (2-21) 

r (^) is the directional reflectance which in this case is the same as the 
instantaneous albedo. 
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Is the bidirectional reflectance. 

Thus, the directional reflectance Is the ratio of the total radiant power per 
unit area W^(0 which is reflected by dA onto a 2it steradians solid angle to 
the total solar radiant power per unit area H(^) incident on dA at a zenith 
angle The bidirectional reflectance involves two directions, namely, 
the direction of the incident radiation which is given by the zenith angle 
C, and the direction of the observer. The latter direction is defined by 
the zenith angle 6 and the azimuthal angle tp. The bidirectional reflectance 
is the ratio of that amount of radiant power per unit area per unit solid 
angle which is reflected by dA onto an element of solid angle dH (6,4*) 
to the total solar radiant power per unit area H(^} incident on dA at the 
zenith angle C . 

The integral over the angle 6 appearing in equation (2-20) is similar 
to the one appearing in equations (2-8) and can be treated in a similar manner. 
That is, one can define I(p) as 


I(p)= p(9»4'»5) sin 6 cos 6d6 

0 


(2-23) 


Then, equation (2-20) can be rewritten as 


A(0 = r(0 = I(p) drD 


(2-24) 


A method proposed for utilizing these equations will be discussed at the 


end of the section entitled Instantaneous/Inversion Technique. 



Instantaneous Net Flux 


Once the representative values of H, W^, and are determined for a 
particular region from a given:' set of n observations 1, taken during a time 
Interval centered at t j , one can determine the net radiant power, or net flux 
Qi(tj) for this region. Hence, Q^Ctj) represents the instantaneous value 
of the net radiant flux for the specific time of day t^ . 


Average Net Flux 


It follows from the above discussion of the Instantaneous energy budget 
that a weighted average (using the configuration factor as a weight) for a 
given region can be' obtained from the m values of computed for m 

sets of observations. 

If, for instance, one is Interested in obtaining the average of Q(^j) 
during a month for a unique value of t^ (say 3:00 PM) from m values, one has 


Q(tj) 


f Fl Ql (tj) 


m 



(2-25) 


Other types of averages which are more suitable for the user’s needs 
can be computed in a similar manner. 

NOTE: In this report only the LWR component will be treated. The SWR 

component will be discussed in a subsequent report. 
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Radiative Equilibrium 


When the temperature of the sensor reaches a steady state, it is said 
that the sensor is in radiative equilibrium. In this state, the power input 
is equal to the power output 0*^ the basis of the definition of 

N(6) for the LWR presented in previous sections, one can write the expressions 
for for both satellites, considered as blackbodles, as 

a 

f ™ NC0(a)3 sin a da (2-26) 

0 




SPHERE P. = A d\p 

in So 


PLATE 


P. = 
in 


A / 


2tt 


dip / ^ NC0(a)Jsin a cos ada 
0 


(2-27) 


where a is the maximum value of the nadir angle, 
m 

The expression for both satellites follows from Stefan-Boltzmann 

law. 


SPHERE 


out 


= S a T 


4 A a T^ 
s 


(2-28) 


PLATE P^^^ - Ag a T (2-29) 

Where, a » 5.6697x10”® watts/ (m^-°K^) ; Stefan-Boltzmann constant. 

It will be assumed only for convenience that the characteristic areas of both 

2 

radiometers is A » Im . Also, N[6(a)] will be written in terms of and 
f [0 (a) ] perieqilati6h:' i('2-7) . 

Then, the above four equations can be rewritten as. 
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(2-30) 


SPHERE 


a 

P. = 2^ fC9(ci)I] sin ada 

in „ 


PLATE 


a 

P. = 2tt ™ fC0(a)I] sin a cos ada 

in 


(2-31) 


SPHERE 


P ^ = 4 a T 
out 


(2-32) 


PLATE 


P ^ = o T 
out 


(2-33) 


Substituting in (2-30) and (2-31) the value of N* as given by (2-9) , one 
obtains , 


/ ^ fC0(oi)3 sin a da 


SPHERE 


p = X! 

in -ir/ 2 


W 


f f(9) sin 0 cos 9 d0 
o 


(2-34) 


PLATE 


/ fC9(a)3 sin a cos a da 

P. = W 

f(0) sin 0 cos 0 d0 ^ 
o 


(2-35) 


But comparing equation (2-5) , that is P » FW, with equations (2-34) and (2-35) 
one sees that the expressions for the shape factors for both radiometers are 
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(2-36) 



a 

• “-f-E© (a-)'3“ sin-a'-da—-* 

f(e) sin 0 cose de 
o 


a 

/ fE6(ct)D sin a cos a da 
s ja 

pi y.ir/2 0 j,Qg 0 ^0 


or, using the definitin of 1(f) given by (2-10), one can write 


(2-37) 


F = 
sp 


_1 . m 

1(f) 0 


fC0(o)D sin a da 


(2-38) 


F t T c\ ^ fC0(a)3 sin a cos a da 

pi 1(f) o 

For the special case of Isotropic radiation 

f [e(a)] " 1 and 1(f) - Then, 


F 


iso 

sp 


2(1 - cos a ) 
in 


(2-39) 


(2--40) 


F 


iso 

pi 


. 2 

sin a 

m 


(2-Al) 


2-13 



TUSDAMENTALS- OF- THE TECHNIQUES 

The only piece of Information that a single satellite observation can 
yield is the total radiant power Intercepted by the radiometer, as given by 
either equation (2-30) for a sphere, or (2-31) for a plate. This power 
represents the siim of those radiances originating at each of the area 
elements within the FOV and directed toward the radiometer. Therefore, it 
is impossible to obtain from a single, wide-angle measurement, and without 
additional information of any kind, an exact or rigorous solution for either 
W or W for any area within the FOV of the radiometer. 

A similar situation exists for measurements by narrow-angle radiometers. 
The only Information that is contained in a single observation of an area 
element is the radiance in the direction of the radiometer on the satellite. 
Scanning by this radiometer only adds information about adjacent area 
elements and nothing additional about the original area element under consid- 
eration. Hence, the fixndamental difficulty one faces when trying to determine 

the W field from a measurement of the radiance N (or W from N ) is that the 
e r r 

total angular distribution of N (or N^) can not be deduced frcao the 

measurement", and hence "'can not be deduced, from the measurement 

either.’ . 

" It. follows from '.the above discussion that it is. Impossible to determine 

the (or W^) field of an area of., any .size - i.e.,. of an- area smaller than the 
FOV, equal to the FOV, or larger than the FOV of the radiometer - unless the 
angular distribution of N is'.'. also made availalbe. On the other hand, it 
is possible that by utilizing representative values of the angular distributions 
of N and obtained from previous observations the problem can be rendered 
solvable in a satisfactory manner. Furthermore, these angular distribution 
values can be refined in the future by Increasing the number of measurements 
of these distributions for different regions of the E-A system, as well as by 
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Improving the accuracy and precision :.of these measurements. 


Therefore, one can utilize an empirical model which portrays the angular 
emitting and reflecting characteristics of the E-A system, based on previous 
observations (satellite, aircraft, etc.) In order to solve the problem. At 
Drexel University, this approach was followed by making certain assumptions 
which are justifiable on the basis of physical processes and observational 
data accumulated by other Investigators (reference 4) . These assumptions are 
discussed In the following subsection. 

Basic Assumptions 

The following basic assumptions were sufficient to develop a technique 
for computing and from radiant power measurements taken by spherical 
and horizontal flat circular radiometers. 

Division Into regions. The surface area of the earth - atmosphere system Is 
assumed divisible Into regions which have emitting and reflecting characteristics 
significantly distinct from those of adjacent regions. The criterion for 
Implementing the division must ultimately be based on results of previous 
satellite observations. The values of and computed for a given region 
are considered to represent the mean of the meteorological variations taking 
place In the time Interval At(e.g., a month) during which the observations were 
made« 

Position and extent of each region . The position and extent of each region 
Is assumed to be available or can be approximated from data gathered by previous 
Investigators. 

Angular distribution of N and N^ . The angular dependence of N for LWR, and of 
for the reflected solar SWR are available, or approximations can be made, 
from previous observations. 
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Although two different techniques for computing the radiant emlttance 
and radiant reflectance fields have been developed at Drexel, the Initial 
computational steps are common to both techniques. These common steps are 
discussed In detail In the following subsection. 

Mathematical Bases 

In this subsection, the mathematical bases which are common to both techniques 
are treated In detail. At the end of the subsection, the principal differences 
In the two techniques for accomplishing the final goal are listed. 

The problem to be solved Is Illustrated in figure 3-1, which depicts a 
satellite (which for simplification Is assximed to be crossing the equator) with 
three different regions (R^, R^, and R^) falling, partially or totally, within 
the FOV of Its radiometer. The following development Is for a spherical 
radiometer; however, the treatment for horizontal flat sensor Is similar. 

A spherical radiometer of unit cross-sectional area is considered. Using 
P <* FW, as given by equation (2-5) , one can write for the total power contributed 
by those portions of the three regions within the FOV of a radiometer on board 
the satellite at position S In figure 3-1 the following. 


P. 

in 



+ P, 


+ P, 


(3-1) 


or 


P. 

in 


F- W , + F- W - + F, W , 
1 el 2 e2 3 e3 


(3-2) 
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That is, by using (2-30) to represent the power originating at region and 
which is intercepted by the satellite, one can write 


P = F, W , 
1 1 el 


’^1 

N? / ^ dij^ / , CQ (u)!]sin a da 

'f'14 «u ^ 


(3-3) 


Where , 

• radiant power from region intercepted by the;, radiometer, 
*■ zenith radiance of region R^. 

■ lower boxmdary for t|i in region R^^. 

B upper boundary for i{* in region R^. 

“ lower boundary for a in region R^^. 

■ upper boundary for o in region R^^. 


To facilitate the discussion, an isotropic radiation field will be assumed. 
Hence, f [6(a)] ■ 1, and one obtains 


P =F. W , 
1 1 el 


„1S0 , lu J , - lu , , 

N- / dip J sin a da 

^ “1£ 


(3-4) 


Where 


- Isotropic radiance of region R^^. 

Using expressions similar to (3-4) for P2 and we can rewrite (3-1) for 
the power input to the radiometer, as 


P=Nf®° / d\f) / sin a da + / dij; / sin a da + J" d\p / sin a da (3-5) 

1 R^ 2 r 2 R2 ^ R3 ^3 
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Where R^, R2> and R^ Imply that the Integration limits are those corresponding 
to the regions 1,2, and 3, respectively. 

Equation (3-5) can be rewritten as 

3 

iso 

P = Z N. // sin a da d li; c\ 

i=l 1 (3-6) 

By using equation (2-12) , P can be expressed in terms of W^^ rather than 
that is, 

P = — Z W . // sin a da d if) (3-7) 

ir i=l ex 

The main problem at this point is to decide which is the most advantageous and 
efficient manner to perform the integration indicated in equations (3-6) or 
(3-7). 

Two main approaches were considered; these are discussed in the following 
paragraphs together with the advantages that each presents. 

The following different methods were considered. 

(a) The Integrations indicated in equation (3-6) could be performed 
numerically (e.g., by using the trapezoidal rule) between the 
boundaries given for each region in terms of longitude and latitude. 

(b) The integrations could be performed geometrically. That is, by 
dividing the surface area of the E-A system Into a large number L of 
area elements, one could simply sum up the radiances that each of 
the area elements emits in the direction of the satellite in order 
to accomplish the integration indicated by (3-6) . 

It was decided that the second approach was more advantageous for the following 
reasons : 
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(a) In the first approach, since the boundaries of the regions do not 
necessarily follow longitudinal and latitudinal lines of constant 


value, the number of separate integrations to be performed could 
often be impractically highi Furthermore, the limits of 
integration for the angles a and if> could Involve difficult tri- 
gonometric expressions which could make the implementation of 
this procedure highly impractical. 

(b) In the second method, once the centroids and boundaries of each of 
the area elements have been accurately defined, the regions can be 
easily defined in terms of a group of adjacent area elements. 

(c) In the second method, a procedure for testing whether a particular 
area element AA should be considered to be within the FOV or not 
can be easily Implemented if the longitude and latitude of the 
centroid of AA is known. If the centroid is on, or Inside of the 
perimeter of the FOV, AA is considered to be within the FOV; 
otherwise, it is considered to be outside the FOV. On the average, 
it is expected that as many area elements will be accepted as will 
be rejected in each integration. Furthermore, since it is at the 
limb that the area elements are tested to be accepted or rejected, 
and since these area elements have much smaller shape factors than 
those situated close to the SSP, it is clear that the errot would 
be negligible even if the number of area elements accepted do not 
match the number rejected. 

Therefore the second method will be fbllowed in which the 
surface of the E-A system is : divided into a large number L of 
area elements AA, and each area- element is identified by the 
longitude and latitude of its centroid. 
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In order to Illustrate the computation of the radiant power Intercepted 
by a satellite radiometer whose FOV comprises more than one region, a FLAT 
E-A system will be assumed. This type of system greatly simplifies the 
results. Only three regions and three FOV's will be considered, which are 
depicted in figure 3-2. The portion of the flat E-A system shown in this 
figure Is comprised between 0° and 70 ° longitude, and between -50° and +50° 
latitude, which Is adequate for the purposes of this discussion. An area 
element Is considered within the FOV If Its centroid lies on, or within, the 
perimeter of the FOV, as per the criterion stated earlier. In figure 3-2, 
the area elements have been numbered In order to Identify them more easily. 

In TABLE 3-1, the elements that are considered to be contained in each of: the 
three FOVAs are tabulated. It Is seen that In this Instance, each FOV 
contains exactly twenty-six area elements. 

The equations used to calculate the radiant power that an area element 
AA contributes to one observation are presented below. The symbols 
entering In these expressions are shown In figure 3-3 which Is a schematic 
of a radiometer at a position- S over the flat E-A system. 

Let A Pyjj the power that the ith area element in the Kth, region 
contributes to the jth observation. ^en,one writes for a spherical radiometer 


AP®., = N.,, AA,., cos 0,,, 
ijk ijk ijk ijk ^2 

ijk 


For a horizontal flat circular radiometer, one writes. 


(3-8) 


AP?,, = N. A A. 
ijk ijk xjk 


cos C( 


cos 6 


ijk 


iik 


ijk 


(3-9) 


N. " Radiance of AAi in the direction of the radiometer 
ijk 


A ~ B Characteristic area of the satellite radiometer, which 

for the sphere and plate is ira^ (a * radius of sphere 
and plate) . 
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TABLE 3-1. 


Observation 

No. 


1 

1 

1 

1 

1 

1 

1 

1 


2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


3 

3 

3 

3 

3 

3 

3 

3 

3 


Data from three satellite observations. 


Region Area Element No. of Elements 

No . No . Elements Per Obs. 


1 201-203 

1 186-190 

1 172-176 

3 158 

2 159-162 

3 144 

2 145-148 

2 131-133 


3 
5 
5 
1 

4 
1 
4 

3 26 


2 159-161 

3 144 

2 145-148 

3 130 

2 131-134 

3 116 

2 117-120 

3 102 

2 103-106 

2 89-91 


3 
1 

1 

4 
1 
4 
1 
4 

3 26 


2 117-119 

3 102 

2 103-106 

3 88 

2 89-92 

3 74 

2 75-78 

3 60-64 

3 47-49 


3 
1 

4 
1 
4 
1 

4 

5 

3 26 
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Figure 3-3. 


Radiance emitted by of * flat E-A system and 

Intercepted by the radiometer S. 
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AA 


ijk 


ijk 


1th area element in the kth region within the 
FOV of the jth observation. 


Zenith angle of the radiometer at AA. 


ijk 


a... ■ Nadir angle of AA. ,, at the radiometer 

IjK ® Ijk 


d. ■ Distance from AA... to the radiometer 
Ijk Ijk 


Let It be assumed that is a function of the zenith angle 6^^^ 

only (i.e., that It Is Independent of the azimuthal angle Then, one 

ijtt 

can write 


Where, 


■*ijk ■ “Uk '(®yk> 


C3-10) 





radiance In the zenith direction at AA 


ijk 


llmb- 

at- 


AA 


darkening function (LDF) or anisotropic fact6.rv 
ijk 


It will be further assumed that a region K Is characterized by the 
following : 

(a) The value of the LWR radiant emlttance In the kth region has a 
representative average value 

(b) The LWR radiance in the zenith direction has also a representative 
average value in the kth region. 

(c) The LDF has a representative functional form throughout 

the i^h region. 

Hence, for the kth region one rewrites (3-10)as 


^ijk “ \ \ (®ijk) 


(3-10’) 
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To avoid carrying A in the following expressions , it is assumed that 

s 

2 

A ■ 1 m . This assumption does not affect the validity of the results, 
s 

2 

By substituting (3-10) into (3-8) and (3-9) , and making A ■ Im , one 

s 


^^ijk“ \ ^k ^®ijk^^ijk ®ijk 


COS a,.- 
Ilk 




In the section BASIC IDEAS, the following expressions were Introduced 

2 

to show the relationship between Nj^ and From equations (2-10) and 

(2-11) one then writes, respectively. 


I(f, ) = f ' f, (9)sin 9 cos 9 d9 
jC 0 


^ y 

k 2 tt I (f^) ek 


Substituting (3-14) into (3-11) and (3-12) one obtains. 


• M ) A A COS 9 . .. 

A p® . T ^ w 

ijk *- 21 (f ) it-* ,2 ek 

ijk 


cos 9 , cos a, ,, 


>P _ fk--ijk- .^AA -— ijk —ijk 
•ijk ^ 21 (f,.) ,2 




ek «-16) 
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Although the subscripts in identify the position, the observation, and 

the region, respectively, of the area element under consideration, these sub- 
scripts can be dropped here since they appear in other variables in the same 
expression. Also, the magnitude of AA has been assumed to be the same for all 

area elements; and hence, it is unnecessary to Identify it with the subscripts. 

2 

It was shown in expression (2-5) that if A « Im , the power P is related 

s 

to the radiant emlttance W by means of the shape factor F, that is 

e 

P = F W (3-17) 

e 

Hence, using the definition of the shape factor, one obtains from (3-15) and 
(3-16), 


pS =r ^ ] 1 °^ 

xjk 


(3-18) 



1.1 - ip A A -] 
21 (f^) 7T ^ 


COS 9 , cos a. ,, 
ijk 1,1 k 


(3-19) 


From figure 3-3, it is seen that for a flat E-A system (which is the type of 

2 

system being used in this development) d^^^ is given by 





(3-20) 


which is different from the equivalent expression for a spherical E-A system, 
which from figure'. 2-1 can be seen to be 


r^^^ = + (R+H)^-2R(R+H)cos y 


(3-21) 
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where , 


""ijk 


R 

r 


distance from 
E-A system. 


radius of the 


angle between 
center of the 


AA. .. to the radiometer in the spherical 


E-A system, 

and the SSP, measured at the 



In order to obtain the power contributed by the kth region to the jth 
observation of both, the sphere and the plate, equations (3-15) and (3-16) 
are summed up over 1. 


p! = s ap! = i w u 

jk i=i ijk i=i ijk ek 


P?, = Z AP?., = Z F^.i, W . 
jk i=i ijk i=i ijk ek 


where , 


(3-22) 


(3-23) 


I * total number of area elements AA. .. of the kth. 

region that were under the FOV oz^lhie jth obser 
vatlon. 


Similarly, the configuration factor F., represents the contribution by 
the kth region to the jth observation and Is obtained as follows. 



I FT., 
1=1 ijk 


(3-24) 



Z F^,, 
1=1 ijk 


(3-25) 


And for the three contributions from the regions Rj^ (k * 1,2,3) to the jth 
observation one obtains by combining (3-22) with (3-24) and (3-23) with (3-25) 


3-15 



(3-26) 


s 3 s 3 s 

p! = E p!. = 2 F®. W , 

J k=l Jk k=l ek 


P? = Z P?, = Z F?, W , 
J k=l 3^ k=l jk ek 


(3-27) 


Or, writing in detail these expressions, one obtains for the sphere. 


P® = F®. W , + F® W , + F® W _ 
1 11 el 12 e2 13 e3 


p® = F® W + F® W + F® W 
2 ^21 el ^22 e2 ^23 e3 


(3-28) 


P® = F® u + F® W + F® W 
3 31 el 23 e2 33 e3 


And for the plate one has. 


>P = pP w + F^ W + F^ W 
1 ^11 el ^12 e2 13 e3 


P? = 


F?, W , + fL W „ + fL W 


2 21 el 22 e2 23 e3 


(3-29) 


pP = F^ W + F^ W + F^ W 
3 31 el ^32 e2 33 e3 


By denoting the matrix of the coefficients for the sphere and plate, respectively, 
as F® and F^, one can express (3-28) and (3-29) in matrix form as follows. 


F { W } = { P® } 

e 


(3-30) 


pP { W } = { pP } 


(3-31) 
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Where, {W^} is the column matrix of the radiant emittance whose elements 
are the We^ (k - 1,2,3) values of the three regions. {P^} and {P^} are 
column matrices whose elements are the individual power measurements by the 
spherical and flat radiometers, respectively. 

The mathematical expressions presented up to this point are common to 
both of the techniques developed at Drexel. From this point on the two 
techniques proceed by different paths. The fundamentals of these two 
techniques will be broadly discussed in the next paragraphs before proceeding 
to the next sections in which each of the techniques and the results obtained 
by applying them are discussed in detail. 

Instantaneous /inversion technique . In this technique, the number m of ob- 
servations matches the number of regions n (l.e., the n unknown values of 
W^) , and hence, the solutions are unique. It is termed instantaneous since 
the procedure is Implemented by taking all the cF'n observations during one 
single pass or orbit of the satellite. Several sets of instantaneous results 
obtained at different time intervals can be grouped together and weight- 
averaged over the complete time period comprising the total number of obser- 
vations. The weights to be used are the configuration factors of those 
region segments appearing within each of the FOV's, or the areas of the 
region segments. 

Therefore, in this procedure, the solution to the n simultaneous 
equations of the type shown in the set of equations (3-22) for n=3 can be 
accomplished by inverting the configuration factor matrix. The inverted 
matrix then, when multiplied by the column matrix of the n power measure- 
ments P. yields the n Wej^ values one seeks. 
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Best-fit/inverslon technique. In this technique, the number of observa- 
tions m is larger than the number of regions n (i.e., the n unknown values 
of W^) . The n simultaneous equations required to solve the problem are 
obtained by using an extension of the method of least squares. These n 
simultaneous equations are solved by using a matrix Inversion subroutine. 

The n We^ values determined in this manner represent the partial coeffici- 
ents of a surface which best fits the data resulting from the observations. 
These represent a mean atmospheric situation portraying the overall condi- 
tion for a significantly large time period (e.g., a month), and no instan- 
taneous results are ever obtained from application of this technique. Here 
again, each of the m observations will be represented by an expression simi- 
lar to those shown in the set of equations (3-26) or (3-27). 

The details of this procedure are presented by means of an illustration 
in the section entitled "Best Fit/Inversion Techniques." 
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INSTANTANEOUS/INVERSION TECHNIQUE 

A mathematical technique for computing the value of will. first be 
presented. Afterwards » the technique will be tested by applying It to a 
simplified case of a flat E-A system. 

Mathematical Development 

The set of simultaneous equations (3-28) for the case of n unknowns can 
be written in matrix form as follows 


{P} » F {W } 
e 


(A-i) 


where 


{P} > column matrix whose elements are the n power 
measurements 


{W } 
e 


column matrix whose elements are the hypothetical 
values of for each of the n regions 


n by n configuration factor matrix whose elements 

are given by (3-18), and which, can be written as, 


F = 


F 

11 

F 

12 ••• 

^in 

F 

21 

F 

A 22 • • • 

^2n 

nl 

^n2 ••• 

F 

nn 


(4-2) 


F the Inverse of F, ,1s then computed in order to solve equations 
(4-1) for by operating with F ^ on {P}, that Is, 


f“^ {P} = {w } 
e 


(4-3) 
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So far, It has been tacitly assumed that F and {P} are exact (l.e., do not 


contain errors). Now, however, it will be assumed that the observations {?} 
include uncertainties {6P}. 

Then, the actual equation to be solved is not the exact equation (4-1) 
but the perturbed equation (reference 6) 


F {W} = {P + 6P} (4-4) 

e 

where {W^} has been introduced to represent the exact solution of the 
perturbed equation (4-4) as opposed to {W^} which is the exact solution of 
the exact equation (4-1) . One can then write 

{w'}= {W + 6W } (4-5) 

e e e 

Again, using the inverse of F one obtains from (4-4) and (4-5), 


-1 


{P + 6P} = {W + 
e 


6W } 
e 


(4-6) 


or. 


f“^{P} + f“^{6P} = {W } + {6W } 

e e 

Subtracting (4-3) from (4-7) one obtains 

F“^{6P} = {6W } 
e 

If the elements of the error matrix {6W^} are large for small error elements 
{6P}, the perturbed equation (4-4) is known as ILL-CONDITIONED, the INVERSE 
MATRIX F”^ is termed UNSTABLE, and the matrix F is called ILL-CONDITIONED 
(references 6 and 7) . 

If the elements of {6W^} are small, or acceptable according to some 
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prescribed accuracy requirements, then expression (4-8) has solved the problem 

of retrieving the values of and nothing else needs to be done. On the 

other hand, If the resulting elements of .{6W } do not meet the accuracy 

e 

requirements for a given set of small uncertainties {6P}, then the problem 
has just begun. This Is the situation that will be treated In the following 
paragraphs . 

Let It he assumed that the elements of {6W } are unacceptable, but that 

e 

the unstable matrix F ^ can somehow be modified to become a stable matrix 
F which, when operating on {6P}, yields a set of new errors in which 
are now acceptable. These errors are the elements of the new column 
matrix {dW^} resulting from the following equations. 

f“^{P + 6P} = {W + 5W } (4^9) 

e e 

or, 

F“^{P} + f“^{5P} = CW } + {6W } (4-10) 

e e 

But it is known that {W^} Is related to the exact power column matrix {P} 
by (4-3) . That is 


{W } = f“^ {P} (4-11) 

e 

Subtracting (4-11) from (4-10) one obtains, 

F"^{P} - F"^{P} + F“^{6P} + {5W } (4-12) 

e 

or 

[f“^ - F”^] {P} + f"^ {6P} = {6W^} (4-13) 
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Equation (4~13) shows that the new Inverse matrix can contribute In the 
following two ways to the value of {6W^}. 

(a) The n by n matrix resulting from subtracting the old Inverse ma- 
trtx from the new one appears ; multiplying the powerimatrlx {P}. 

(b) The new n by n Inverse matrix F ^ multiplies the power uncer- 
tainty ; matrix {6P}. 

In the subsection entitled "Matrix Stabilization,'", a detailed account of a 
scheme found to stabilize the Inverse matrix F ^ Is given. In that sub- 
section Is also explained the physical significance of modifying the unsta- 
ble Inverse matrix. 

In order to Implement this technique, a computer program would be 
required for performing the numerical Integrations, matrix Inversions, matrix 
stabilizations, and tests that would be deemed necessary. Furthermore, 
another computer program would be required for dividing the surface of the 
spherical earth-atmosphere system Into a finite number of area elements. 

The output data of the latter program would be used In the former for 
accomplishing the numerical Integrations. 

Nevertheless, before engaging In writing these large and sophisticated 
computer programs. It was decided to subject the technique to a test. 

By considering a hypothetical flat earth-atmosphere system, some of the 
equations, as well as the overall computer program could be greatly 
simplified. On the other hand, the procedure developed to stabilize the 
Inverses of the matrices would be rigorously tested. In the next subsection, 
the flat E-A system is discussed. 
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Flat Earth Application 


By considering the E-A system to be flat rather than spherical, the task 

vas greatly simplified and less time consuming, x^i^'e 'the results obtained 

can be Interpreted as a representation of actual physical situations In a 

spherical system. This flat E-A system was assumed to be a rectangle of 

360*^ of longitude by 180° of latltutde. The longitude Is measured westward 

(to the left) from 0°, at a hypothetical Greenwich meridian, to 360°. The 

latitude is measured northward (upward) from 0°, at a hypothetical equator, 

to 90°; and southward (downward) from 0° to -90°. It was found convenient, 

for interpreting some of the results, to consider the flat E-A system to 

be tangential to the spherical E-A system at the SSP, as shown In figure 

4-1. At any point in this flat E-A system, one degree of longitude is" the same 

as one degree of latitude and each Is equal to 100 km. The surface area 

of this system was divided into a total of 2592 equal area elements AA of 

5° by 5°. The area of an area element Is given by AA » 500 km x 500 km » 

2 11 2 

250,000 km « 2.5x10 m . The surface area of the system was also divided 
into 162 regions of 20° by 20°, or 2000 km by 2000 km. It is to be noted 
that the area elements are used to perform the numerical or geometrical 
summations or integrations. Each area element AA Is Identified by the 
longitude and latitude of its centroid, and hypothetical values of and 

albedo can be assigned to AA. However, In the present report only Is 

included. Direct and reflected SWR will be discussed in a subsequent report. 

A region is characterized by a representative value of (and/or W^(0)), 
and variations of (or W^(0)) within it are neglibibly small or undetectable 
by a radiometer. 

The radius p of that area of the flat E-A system within the FOV of the 
radiometer can be quickly claculated with the aid of figure 4-1*. From this 
figure. It can be seen that 
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Figure 4-1. Flat E-A system parameters. 
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(4-14) 


and. 


Hence, 


a 

m 


sxn L 


R + H 


p = H tan a 
m m 


= H tan Csin ^ ^ )3 

For a height H ■ 800 km, a » 62.74°, and 
” m 


(4-15) 


(4-16) 


p (H=800 km) = 1550 km (4-17) 

m 

Since in the flat E-A system 1° « 100 km, one has 

p’ (H=800 km) = 15.5° (4-18) 

m 

As shown in figure 4-1, p is the distance between the SSP and the centroid 

of the area element AA under consideration. Therefore, according to the 

criterion previously discussed for deciding if a given AA^ should be 

considered within the FOV or not, one should compare the corresponding value 

p. with p . If p < p , AA, is considered to be within the FOV. Therefore, 
1 la 1 TO 1 

one can write into the computer program a simple scanning scheme for testing 
all area elements within the area of interest. Of course, the total of 2592 
AA's can be scanned each time a particular AA is considered for Inclusion 
in the computations, but computer time is saved if one can program the 
boundaries of an area somewhat larger than the area of Interest. This point 
Is Illustrated in detail in the following exercise which was actually used 
to evaluate the validity of the technique being discussed. 
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Error F.ree Observations 


In this subsection » a spherical radiometer Is simulated in a straight line 
trajectory 800 km above the E-A system. A horizontal flat circular radiometer 
will be asstjmed to be coincident with the spherical radiometer at all times 
and hence the surface arear;of the E-A system intercepted by both F0V*s will 
be identical. Arbitrarily, six regions and six satellite positions have 
been selected which are portrayed in figure 4-2. Although a total of fifteen 
regions are shown in this figure, only six of then are observed at one time 
or another by the satellites; and hence, only the values for these six 
regions needed to be shown. Nevertheless, the values of eight regions 
were actually included. The six satellite positions are identified by num- 
bers one through six as the satellites travel from south to north in their 
common trajectory. The longitude and latitude of the SSP for each observation 
is shown in parentheses next to the number of the observation. The longitudes 
and latitudes of the centroids of all area elements can be easily found by 
referring to the longitudes and latitudes indicated at the margins. For 
example, the longitudes and latitudes of the centroids of the two area ele- 
ments positioned north and south of SSP N 0.4 are as follows. 

North; long. ,22.5**, lat. +12.5**. South: long. 22.5**, lat. +7.5**. All 

the data pertaining to those area elements contained within the eight regions 
bounded by longitudes 0** and 40**, and by latitudes -20** and 60** were fed into 
the computer program generated to test the technique. 

In order to illustrate how the results for the situation portrayed in 
figure 4-2 . were obtained, one may consider the area element AA(12.5,32.5) , 
that is the area element whose centroid has longitude 12.5** and latitude 
32.5**. The perimeter of this area element is shown marked with a broken line 
in figure 4-2.. For purposes of this illustration, satellite position No. 6 
which has a SSP at 24** longitude and 20** latitude is selected. The first 
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Figure 4-2* Six Satellite FOV’s over six regions of the flat E-A system. 
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step Is to compute (12.5,32.5) for this situation and compare it with 

to determine if ^A (12.5,32.5) is or is not within the FOV of the satellites 

at position No. 6. In this notation the pg (12.5, 32.5) means the distance 

(in m) between SSP No. 6 and the centroid of M (12.5,32.5); while p. means 

oin 

the radius of the FOV for SSP No. 6, which is the maximum value p can attain 
for observation No. 6. 

As shown in the previous subsection for H » 800 km, p^^ » 1.55xl0^m. 

Using primes to Indicate the distances in degrees . one calculates p^ (12.5,32.5) 
as follows. 


where. 



II 

Q. 

long 

) * 
lat'' 

A? 

long 

- (12.5-24.)^ = 

(-11.5)^ 

- 132.25 


- (32.5-20)^ = 

(12.5)^ 

- 156.25 

P* 

- distance from 

SSP to AA in degrees 


(4-19). 


Then, 

p’ - (288.5)^ 17.0“; = 15.5° (4_2o) 

f 

Hence, p’ > p^ .which means -that the area element A A (12.5, 32.5) is NOT 
o bm 

within the FOV of the radiometer at position No. 6. 

In order to calculate the shape factor and the power Increment contributed 
by an area element (which must lie within the FOV of the satellites) an area 
element lying inside the FOV of observation No.l is selected. This element 
is AA (7. 5, -12. 5) and its boundaries are also marked with a broken line in 
figure 4-2. Since the whole trajectory is at a height of 800 km, the value 
of for position No.l is the same as Nbi6, i.e., p^^ = 15.5°., The calculation . 
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(4-21) 


of (7. 5, -12. 5) is then, 

aJ = (7.5 - 19)^ = 132.25 
long 

A^^^ = (-12.5-(^5))^ = 56.25 
.= (188.5)^ = 13.7° 

Thus, p'j^< p |^ and AA (7.5, - 12.5) contributes radiant power to the radiometers 
of the satellites at position No.l. 

In order to simplify the computations of this run, the radiation field 
was assumed isotropic. This means that in equations (3-18) and (3-19), 
f(ep » 1 and 1(f) » 

Hence, for a spherical satellite the shape factor for the ith area element 
in the k'th region and within the FOV of the jth observation is 


. . COS0 . 

„s _ AA ijk 

^ijk - — ^2 

‘^ijk 

From figure 4-1, one obtains for the flat E-A system, 

o H 

COS 0 


"Uk 

,3 ^ 2 .3/2 

‘^ijk = (H + Pijk> 


(4-22) 


(4-23) 


Hence, 


AA 


H 


^ijk " ir (HZ + 


For the case \mder consideration p. ,, *» 13.7x10 m, and 

ijk 


AA » 2.5xl0^^m^ 


H - 8x10 m 


Then, 


(4-24) 
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By (3-24) all the shape factors of the area elements within region No.l, and 
within the FOV of SSP No. 1 are added to yield the first element of the 
configuration factor matrix. The value of this element was calculated by the 
computer program TARA 7-15-1, and Its value, shown in figure 4-3, Is 

F®i = 0.484847428 (4-25) 

Then, according to (3-17) for a spherical radiometer with a characteristic 

2 21 s 

area A^ va 1 m , one can write for the power ^11 which region No.l 

contributes to observation No.l 



(4-26) 


(4-27) 


In a similar manner, the radiant power contributed by regions 2,3, and 4 
are computed and added to get the total radiant power Intercepted by the sphere 
at position No.l. It Is noted that since regions 5 and 6 do not appear within 
the FOV, their contributions are zero. Then, one writes. 


P =AP +AP +AP +AP 
^1 ^11 12 13 14 


(4-28) 


P® = F,,W , + F,_W - + F-«W _ + F^,W , 
1 11 el 12 e2 13 e3 14 e4 


(4-29) 


This power was calculated by the computer program and its value, shown in 
figure 4-3, is 


P® = 262.892068914 w (4-30) 
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An identical procedure would be followed In order to obtain for the 
horizontal flat circular satellite. There Is only one difference between the 
two calculations which can be seen from comparing the expressions for the 
shape factor of the sphere (S-IS) and for the plate (3-19) . The latter has 
an additional factor, namely, cos angle Is the nadir 2 mgle of 

AA^^^ measured at the satellite, and cos is used In order to accouht': 

for the projection of the plate area onto a plane perpendicular to the direction 
of the radiance N. (refer to figure 3-3). In the flat E-A system. 

IjK 

and hence, cos \jk* (3-19), the expression 

for the shape factor for the plate, assuming an isotropic radiation field. Is 


ijk 


2 o 

. . cos 9 . 

A A Ilk 

j2 

•"idk 


(4-31) 


The values of the variables In this expression are the same as for the 
sphere. Using the value of cos given by (4-23) then, one has. 


.P = 
Ijk 


A A H 
ir ,4 


(4-32) 


ijk 


F?., = 0.008 


(4-33) 


In this case this value is one half that of the sphere; however, this ratio 
varies for different area elements as would be expected. The area element 
here considered Is at the limb where cos attains Its smallest value. For 

an area element at the nadir position, cos = 1, and the shape factors for 

both sphere and plate have Identical values. 

Again here, adding all the contributions from region No.l to the plate In 
position No.l, one obtains the value of Fj^^ shown in figure 4-4, that is. 
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F^i * 0.378000655 


(4-34) 


Then,; again, according to (3-17) for horizontal flat circular radiometer with 

2 2 

a characteristic area A “ ra - Im^ one can write for the power P,? which 

8 11 

2 

region No.l (W ^ • 236.0 w/m ) contributes to observation No.l 


= 89.20815458 w ^^“35 

Again here, the powers contributed by regions 2,3, and 4 are calculated 
in a similar manner and added in order to obtain the total radiant power 
intercepted by the plate at position No.l. Since the FOV's of the sphere 
and plate coincide, regions 5 and 6 must not appear in the FOV of the plate 
at position No.l either. The expression for P^ equivalent to (4-29) is 



F W 
Ik ek 


(4-36) 


The value of P^ calculated by the computer program TARA 7-15-1 is shown 
in figure 4-4; this value is 


pj » 190.073940996 W 

In the previous paragraphs , it has been shown in detail how to calculate 

the elements of the configuration factor matrices F 6ind power column matrices 

{P} for spherical and horizontal flat circular radiometers. The elements of 

the column matrix {W^} were assumed known in order to simulate and compute 

the radiant power observed by both satellites. In this manner, the elements 

of the column matrices {P} for both radiometers were calculated by the 

computer program. Figures 4-3 and 4-4 show all the elements of the square 

F matrix and of the two column matrices {W } and {P}'. These figures show 

e 

theimatrix. products ;fofr the sphere and plate, which can be written symbolically 
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C0FSPH(K,1) 

C0FSPH(K,2) 

C0FSPH(K,3) 

C0FSPH(K,4) 

0.484847428 

0.384856899 

0.127704158 

0.108360067 

0.276829292 

0.276829292 

0.276829292 

0.276829292 

0. 108860067 

0.127704158 

0.384856899 

0.484847428 

0.017068977 

0.035500329 

0. 355753027 

0.595573384 

0.000000000 

0.000000000 

0.247608706 

0.596414856 

0.000000000 

0.000000000 

0.137012641 

0.415871635 


C0FSPH(K,5) 

G0FSPH(K,6) 

WE 

PSPHER(K); 

t .a 

0.000000000 

0.000000000 


236.0 


262.892068914 

0.000000000 

0.000000000 


238.0 


264.648803387 

0.000000000 

0.000000000 


240.0 


265.665298558 

0.017068977 

0.035500329 


242.0 


254. 884753820 

0.068729515 

0.143711947 


244.0 


255.881625193 

0.137012641 

0.415871635 


246.0 


269.259475929 

) 


Figure 4-'3. Matrix product computed for the sphere 
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GOFHOR(Ka) 

C0FH0R(K,2) 

0.378000855 

0.279486438 

0.200205230 

0.200205230 

0.063693275 

0.078892673 

0.008805215 

0.018591209 

0.000000000 

0.000000000 

0.000000000 

0.000000000 


C0FH0R(K ,3) 

G0FH0R(I(,4) 

0.078892673 

0.063693275 

0,200205230 

0.200205230 

0.279486438 

0.378000655 

0.251441001 

0.471341713 

0.172299495 

0.473440034 

0.087326867 

0.312709653 


C0FH0R(K,5) 

COFHOR(K ,6) 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.008805215 

0.018591209 

0.040187252 

1 

0.091648780 

0.087326867 

0.312709653 


We PHORIZ(k) 


236.0 


'90.073940996 

238.0 


L91. 396199412 

240.0 


L93. 360972805 

1 242.0 


187.635182979 

244.0 

1 

188.275656474 

1 246.0 


194.868514638 


Figure 4~4. Matrix product computed for the plate. 



as 


{P} = F {W } (4-37) 

e 

The first test to which the technique was subjected Is the following. 
Consider the error-free power measurements contained in the column matrices 
entitled PSPH£R(k) and PHORIZ(k) In figures 4-3 and 4-4 respectively. 

Compute the six values of by operating with the Inverses of the configu- 
ration factor matrices on the power column matrices as Indicated by equation 
(4-3), that Is, by using 


F"^{P} = {W } (4-38) 

e 

For the case In which the observations are free of uncertainties, the six 
values of retrieved should be the same as the six hypothetical values 
of originally given. 

The Inverses F ^ of the original configuration factor matrices for 
the sphere and plate are presented In figures 4-5 and 4-6, respectively. It 
should be noted that the Inverses appear In the transposed form which Is the 
order In which the computer, ordinarily, prints out matrices. 

The six values obtained when performing the matrix multiplications 
(for the sphere and the plate) given symbolically by (4-38) were identical 
to the six original hypothetical values to at least eight decimal places, 
for both satellites. 

The technique was then subjected to the following test. The observations 
were assumed to Include gausslan random uncertalntltles. The procedure 
followed and the results obtained are described In detail In the next 
subsection. 
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4-18 


.1264066378 02 
-.299509655E 02 
.322632341E 02 
-.195404147E 02 
.491022939E 01 
-.287768318E-01 


-5.136852019E 02 
.403290804E 02‘ 
-.463025861E 02 
.286083638E 02 
-.718887652E 01 
.421310440E-01 


,1285398348 01 
-.137579890E 02 
.339027825E 02 
-.286083638E 02 
.718887652E 01 
-.421310441E-01 


-,240860110E 00 
.699220887E 01 
-.198634305E 02 
•195404147E 02 
-.491022930E 01 
.287768319E-01 


T.093927101E 01 
-.191839507E 02 
.103020142E 03 
-.145767511E 03 
.833969093E 02 
-.163760570E 02 


,210357918E 01 
.386080888E 01 
-.252470553E 02 
•379092813E 02 
-.249340703E 02 
.778493063E 01 


Figure 4-5. Inverse of original configuration factor matrix, in ttatvsposed form, for the sphere. 



.131424995E 02 
-»292298882E 02 
.323672225E 02 
-.20X552571E 02 
.499057827E 01 
-.264366054E 00 


-.146981203E 02 
.422074149E 02 
-.498062519E 02 
.315809718E 02 
T.781966270E 01 
.414231230E 00 


.269545375E 01 
-.157270628E 02 
.378035854E 02 
-.315809718E 02 
.781966270E 01 
-.414231229E 00 


-.113983297E 01 
.774441065E 01 
-.203645559E 02 
.201552571E 02 
-.499057827E 01 
,264386054E 00 


.272309581E 01 
-.445058032E 02 
.152739147E 03 
-.209818167E 03 
.120475873E 03 
-,228348995E 02 


-.373343722E 00 
.907613661E 01 
-.328461452E 02 
.472575341E 02 
-.308370566E 02 
.942600877E 01 


Figure 4’-6, Inverse of the original configuration factor matrix, in transposed form, for the plate^ 



Perturbed Observations 


In order to generate the gausslan random errors needed to perturb the 

six radiant power measurements, a computer program using two subroutines 

(BANDUr and GAUSS) was generated. The sigma value of the gausslan 

2 

distribution Is o « 0.5 W/m . The gausslan errors e were added to the 

g 

power values originally computed (l.e., the unperturbed powers) in order 
to obtain the gausslan-perturbed power matrix that is. 


{P} + {e } = {P + e } = {P } (^“35) 

g g g 


where 


P^ > Gausslan perturbed power elements of the new matrix. 


The test to which the technique will now be subjected consists In carrying 

out the multiplication of the configuration factor matrix by the new column 

power matrix {P }. This multiplication Is similar to that in (4-38), except 
g 

that in the latter, the original power column matrix {P} Is used Instead of 

{P }. Hence, the old elements of {W } in (4-38) will be modified by 
g ® 

uncertaintltles c resulting from the uncertalntltles e in the elements of 
wg ® g 

{P }. That is. In matrix form, 
g 


{W } + {e } = {W + e } = {W } 
e wg e wg eg 


(4-40) 


and the matrix product to be calculated can be written as 


{P } = {W } (4-41) 

Where are the Inverse matrices given In transposed form for the sphere 
and the plate in figures (4-5) and (4-6), respectively. The first three 
results obtained by adding Gausslan random perturbations to the exact power 
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measurements are listed In TABLE 4-1. There are three sets of data In this 

table. For each set, the gausslan error for each of the regions appears on 

the first row of the set. The second row contains the uncertaintltles e 

wg 

in the values of for the six regions, resulting from the perturbations 

e of the power measurements of the spherical radiometer. Similarly, the 
S 

third row contains the corresponding uncertaintltles ^ in the values of 
for the six regions, for the flat radiometer. Each of the quantities in 
the last column represents the root-mean-square (rms) values of the six 
errors shown in the corresponding row. Two important facts are to be 
noticed from the last column of TABLE 4-1; these are 

(a) The plate consistently exhibits a larger error than the sphere. 

(b) The results of the matrix inversion are completely unacceptable 
when the observations include gausslan uncertainties. 

At this point, based on the results stated in item (b) above, the 

applicability of this technique was questioned and therefore, it was decided 

to investigate the possibility of improving the situation. However, before 

going into a discussion of how the problem was solved, it will prove 

helpful to Introduce now the results obtained when the observations Included 

systematic uncertainties only, as well as when combinations of systematic 

and gausslan uncertainties were Included. 

TABLE 4-2 shows the errors e in W when systematic errors e were 

ws e •' , s 

added to the power elements of {P}. In this table, two sets of data are 

presented. Each set has in the first row, the systematic errors e for the 

s 

six regions. The second and third rows, as before, contain the errors 

in W for the sphere and plate, respectively. Systematic errors e of equal 
0 s 

magnitude and opposite sign produce errors of equal magnitude but opposite 

sign. 

Two conclusions can be made from TABLE 4-2: (a) Systematic errors in 
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TABLE A-1. Results obtained from observations having gausslan uncertalntltles. 


Set No. 

Type of ; 

Error 

Region 1 

Region 2 

Region 3 

Region 4 

Region 5 

Region 6 

rms of All 
Columns 

1 

Gaussian 

OQ 


1.1430 

- 0.3780 


0.0730 

0.7630 

- 0.3030 

0.7480 

0.6693 


sphere , 

e * 

wg 

- 

17.1903 

21.0556 

- 

17.8317 

12.6009 

- 126.0356 

36.5961 

55.4374 


..Plate, 

e 

wg 

- 

18.6985 

23.9852 

- 

21.1519 

13.9772 

- 188.8153 

47.0500 

81.0685 

2 

' Gaussian 

' 8 

- 

0.6290 

- 0.0900 


0.0220 

1.1540 

0.0010 

0.5750 

0.5869 


Sphere, 

e 

wg 

- 

27.1069 

36.9908 

- 

31.8555 

21.6465 

- 169.1907 

45.9727 

75.6367 


Plate, 

e 

wg 

- 

28.3301 

41.0255 

- 

36.1231 

22.9781 

- 249.4868 

58.6197 

108.0112 

3 

Ga'usslan 

» ^ 

g 

- 

0.6060 

- 0.0060 

- 

0.5100 

- 0.1780 

- 0.0350 

0.3750 

0.3654 


Sphere, 

e 

wg 

- 

20.6392 

26.8407 


13.1619 

6.9388 

- 31.3333 

8.6223 

20.0993 


Plate, 

£ 

Wg 

- 

20.9824 

28.8626 

- 

15.6265 

7.7164 

- 54.7122 

13.1255 

28.1151 

* 

e Is In 
g 

watts 



** 

e 

wg 

Is in w/m^ 

• 
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TABLE 4-2. Results obtained from observations having systematic errors 


Set No. 

Type of Error 

Region 1 

Region 2 

Region 3 

Region 4 

Region 5 

Region 6 

riits of All 
Columns 

1 

Systematic, e^* 

- 0,3 

- 0.3 

- 0.3 

- 0.3 

- 0.3 

- 0.3 

0.3 


Sphere, ®ws** 

- 0.0882 

- 0.5409 

0.0094 

- 0.4641 

0.5549 

- 0.4432 

0.4123 


Plate, e 

ws 

- 0.2552 

- 0.5636 

- 0.1789 

- 0.5007 

0.3662 

- 0.5109 

0.4205 

2 

Systematic, 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 


Sphere , e 

ws 

0.2646 

1.6226 

- 0.0283 

1.3922 

- 1.6648 

1.3297 

1.2370 


Plate, e 

ws 

0.7657 

1.6907 

0.5368 

1.5022 

- 1.0987 

J 

1.5328 

1.2614 


2 

e Is In watts ** e Is In w/m 

8 ws 



the observations do NOT produce large errors in W ; and (b) the plate 

e 

shows larger errors than the sphere. 

Combinations of systematic and gausslan errors e were Included In the 

gs 

power measurements In order to compute the errors resulting In the 

values. As expected, these e errors exhibited more sensitivity to the 

wgs 

gausslan than to the systematic uncertainties. The results obtained are 
shown In TABLE 4-3. The results In this table again show that the sphere 
yields better results than the plate when gausslan random uncertainties 
are Included In the power measurements. So far, however, the two radiometers 
yield acceptable results only when the errors are systematic . 

Thus, as can be seen from the data In TABLES 4-1 and 4-3, the gausslan 
errors are highly magnified through the matrix multiplication of the Inverse 
of the configuration factor matrix and the perturbed power column matrix. 

This Instability of the Inverses of both configuration factor matrices 
(sphere and plate) was corrected as explained In the following section. 

Matrix Stabilization 

Prior to describing the scheme developed for stabilizing the two matrices. 
It Is necessary to return to the topic of Ill-conditioned perturbed equations 
in order to Introduce the concept of the CONDITION NUMBER of a matrix. 

Consider the set of n simultaneous equations represented In (4-1) in matrix 
form as 


F{W } = {P} 
e 


(4-42) 


where 
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TABLE 4—3. Results obtained from observations having combinations of gausslan and systematic errors. 


rms of All 

Set_No. Type of jrror Region 1 Region 2 Region 3 Region 4 Region 5 Region 6 Columns 


1 

Gauss f 

Systematic, 

Total 

e * 

e * 
gs 

— 

1.1430 - 
0.3000 
0.8430 - 

0.3780 

0.3000 

0.0780 


Sphere 

e 

Wgs 

- 

17.1021 

21.5965 


Plate 

E 

Wgs 


18.4433 

24.5488 

2 

Gauss , 

Systematic, 

Total 

e 

gs 

- 

1.1430 - 
0.6000 
0.5430 

0,3780 

0.6000 

0.2220 

4-25 

Sphere 

e 

wgs 

- 

17.0139 

22.1374 


Plate 

e 

wgs 


18.1881 

25.1123 

3 

Gauss , 

Systematic 

Total 

E 

8 

C 

gs 


1.1430 - 
0.9000 
0.2430 

0.3780 

0.9000 

0.5220 


Sphere 

£ 

Wgs 

- 

16.9257 

22.6782 


Total 

e 

wgs 


17.9328 

25.6759 

is £ 
^8* 

e , and e are 
8* gs 

in watts 



** e 

wgs 


0.0730 

0.3000 

0.3730 

0.7630 - 
0.3000 
1.0630 - 

0.3030 

0.3000 

0.0030 

0. 7480 
0. 3000 
1.0480 

0.6693 
0.3000 
0. 7170 

17.8411 

13.0649 - 

126.5905 

37.0394 

55.7449 

20.9730 

14.4779 - 

189.1815 

47.5609 

81.2855 


0.0730 

0.7630 - 

0. 3030 

0. 7480 

0.6693 

0.6000 

0.6000 

0.6000 

0.6000 

0.6000 

0.6730 

1.3630 

0.2970 

1.3480 

0.8718 

- 17.8506 

13.5290 - 

127.1454 

37.4826 

56.0537 

- 20.7940- 

14.9786 - 

189.5478 

48.0718 

81.5040 


0.0730 

0.7630 - 

0.3030 

0. 7480 

0.6993 

0.9000 

0.9000 

0.9000 

0.9000 

0.9000 

0.9730 

1.6630 

0.5970 

1.6480 

1.0890 

- 17.8600 

13.9931 - 

127.7003 

37.9258 

56.3639 

- 20.6151 

15.4793 - 

189.9140 

48.5828 

81.7242 


Is In w/m^ 



F « n by n configuration factor matrix 

{W } B column matrix of n W elements 
e e 

{P} B column matrix of n P elements 

It Is possible that these quantities may include errors resulting from 
either one, or both, of the following (reference 6). 

(a) The data (l.:e., the observations {P}) are inexact. 

(b) Rounding errors are generated during computations. 


Samples were rtin with error- free observations and. all the W values 

e 

retrieved were accurate to at least eight decimal places. Consequently, one 
is jiistlfied in assuming that computational errors have not caused diffi- 
culties, and that only errors of type (a) should be treated in this 
discussion. 

The perturbed equation to be solved is of the form given by (4-4), that is, 

F{W'} = {P + 6P} (4-43) 

where, as indicated in (4-5), 


{W*} = {W + 6W } (4-44) 

e e e 

In this expression, the elements of {6P} are considered to be either 
gausslan errors or systematic errors e^, or a combination of both 

'It is assumed for purposes of this section that the perturbed equation 
(4-43) is ill-conditioned, that is, that the inverse matrix F ^ is unstable. 
It can be shown that for ill-conditioning of this type the following rela- 
tionship exists (reference 6) . 


rel {W } = *rel {P} 
e 2 


(4-45) 
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where 


rel{W. } = 1 

IlswJ 

1 /llWgll 

(relative error) 

(4-46) 

rel{P} = 1 

ll«pll 

/IIP II 

(relative error) 

(4-47) 

S = 1 

l|F II 

• I|f"^ 

(condition number 

(4-48) 


of F) 


The symbol | | j ( denotes the column or row norm of the matrix Inside it, as 
will be explained below. 


Condition number. - Even though there are several definitions of the 
condition number of a matrix (references 6 and 7) only two of them will be 
Introduced In this report. 


max A1 


C, = 


1 . .. 
mxn Ai 


(4-49) 


where 


C2 = l|Fl 


max] Ai j 
min] Aij 

II Fll 


,-li 


(4-50) 


largest modulus eigenvalue of F 


smallest modulus eigenvalue of F 


Column or row norm of F defined, respectively, 
as (reference 5) 


If] I = max 2 If^.I 
J- J 1 


(4-51) 


I|fIL = 


max 2 If. . I 
1 3 


(4-52) 


For the Ill-conditioned slttiatlon considered here and represented by the 
perturbed equation (4-44) , It can be seen from (4-46) that 111-conditlonlng 



depends on the size of the condition number and on the relative error 
rel {?}. Thus, for a given set of power errors {5F}, one would expect 
that If the original Ill-conditioned matrices could somehow be transformed 
Into two well-conditioned matrices, the condition numbers of the latter 
would be smaller than those of the former. 

A computer program was used to obtain the eigenvalues of the two 
configuration factor matrices and to compute the values of C'| for both 
matrices. This program also verified that the eigenvalues computed are 
correct by making use of the following relations (reference 8) 

TrF ■ S X, 

1 ^ 

1f| - ir X 
1 

where, 

TrF « Trace of the matrix F 
X^ “ 1th eigenvalue of F 
|f| “ The determinant of F 

The values of the condition numbers C 2 for both of the original Ill- 
conditioned matrices were calculated. These values are shown In TABLE 4-A 
below. 

TABLE 4-4. Condition numbers of the two original matrices. 

Type of Sensor C 2 

Spherical 131.6 693.9 

Horizontal flat circular 126.4 - 684.7 
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After considerable effort, a scheme was found which rendered both configu- 
ration factor matrices well-conditioned and their corresponding values of 
and became much smaller than those shown in TABLE 4-4. 

Stabilization procedure .- Essentially, this prodedure consists in trans- 
lating the smallest elements in the six by six configuration factor matlrx. 
The translation of each element is performed along the row of the element 
to the position of the diagonal element in that row and it is added to the 
diagonal element. In this manner, the sum of the elements in the row in 
question is preserved. This means that the sum of the row elements always 
adds up tb the configuration factor of the total FOV, which is a desirable 
feature, as will be explained below. 

Referring to the two original matrices in figures 4-3 and 4-4, for the 
sphere and the plate, respectively, one sees that the elements F^^^ and F^^ 
are the lowest value elements in both matrices. For the sphere, these 
two elements have the common value 0.017068977; while for the plate the 
common value is 0.008805215. In each of the matrices, these two elements^ 
were translated and added to the diagonal element Che sphere, 

the original value of this element was 0.595573384; while for the plate, 
the value was 0.471341713. The new values of F^^ for the sphere and plate 
become, respectively, 0.629711338 and 0.488952143. The physical meaning" of 
this element translation can be best seen when evaluating the error it 
introduces into the power calculations. The translation of an element means 
that the element will appear multiplying the value of the column 
corresponding to the diagonal element rather than the value corresponding 
to the column where the element originally appeared. As an example of the 
magnitude of the error introduced, the case of element F^^ is calculated. 

W - = 244.0 W/m^, and W . » 242.0 W/m^. Hence, AW = 242.0-244.0= - 2.0 W/m^. 
and the power error AP introduced Is approximately AP = (-2.0) (0.017), or, 

AP = -0.034 W, which is certainly a negligible error. Physically, this 
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error is equivalent to saying that the radiometers at position No. A looked 
a little more at the limb of region No. 4 and did not see region No. 5 at all. 
Figures 4-7 and 4-8 show the new stabilized matrices for the sphere and 
plate, respectively. The captions of these figures Include the title of 
the computer program (TARA 7-15-1) which stabilized the matrices and then 
used them to compute the values of The stabilization scheme used In 

this progreun establishes lower limits for the- magnitudes of the matrix 
elements for each of the t^b matrices any element whose value is below the 
lower limit assigned to Its matrix Is to be translated as previously des- 
cribed. The limits determined after several trials were 0.032 for the 
sphere and 0.016 for the plate. A different computer program (TARA 7-15-2) 
used different lower bounds as the criteria to carry out the translations. 
These limits were 0.04 for the sphere and 0.02 for the plate. The results 
obtained with the latter limits were not as good as those obtained with the 
former. 

The errors in obtained with computer program TARA 7-15-1 by \ 

first using the two original Ill-conditioned matrices and later using the 

new well-conditioned matrices are presented In TABLE 4-5. Three different 

data groups are tabulated there. These groups refer to the first set of data 

appearing in each of TABLES 4-1, 4-2, and 4-3, as is Indicated in the 

first column of TABLE 4---5. These three groups of data correspond to the three 

types of uncertainties previously Introduced, Guasslan e , systematic 

S 

and combinations of gausslan and systematic The type of Inverse 

matrix used (original or stable), as well as the class of radiometer 
(sphere or plate). Is specified for each row of errors. 

The striking differences exhibited by the results compared in TABLE 4-5 
Is indicative of the effect that stabilization of the configuration factor 
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C CFSPH(K ,1) 

C0FSPH(K,2) 

0.484847428 

0.384856899 

0.276829292 

0.276829292 

0.108360067 

0.127704158 

0.017068977 

0.035500329 

0.000000000 

0.000000000 

0,000000000 

0.000000000 


C0FSPH(K;J) 

C0FSPH(K,4) 

0.127704158 

0,108360067 

0.276829292 

0.276829292 

0.384856899 

0,484847428 

0.355753027 

0.595573384 

0.247608706 

0.596414856 

0.137012641 

0.415871635 


CflFSPH(Ki5) 

c:qfsph(k;6) 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.017068977 

0,035500329 

0.068729515 

0,143711947 

0.137012641 

0.415871635 


Figure 4-?7. Stabilized matrix for the sphere. Computer program TARA 7-15-1 
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C0FHOR(K,l) 

C0FH0R(K,2) 

0.378000655 

0.279486438 

0.200205230 

0.200205230 

0.063693275 

0.078892673 

0.000000000 

0.018591209 

0.000000000 

0.000000000 

0.000000000 

0.000000000 


Figure 4-8. Stabilized matrix 


C:0FH0R(K,3) 

C:OFHOR(K,4) 

0.078892673 

0.063693275 

0.200205230 

0.200205230 

0.279486438 

0.378000653 

0.251441001 

0.488952143 

0.172299495 

0.473440034 

0.087326867 

0.312709653 


CCFH0R(K,5) 

COFHOR(K;6) 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.018591209 

0.040187252 

0.091648780 

0.087326867 1 

0.312709653 


for the plate. Computer program TARA 7-15-1. 



TABLE A~5. Comparisons of results of original and stabilized matrices for the first set of data shown in TABLES 
• 3-2, 3-3, and 3-4. 


Table 

Type of Error 











rms of All 

No. 

Type of Matrix 


Region 1 

Region 2 

Region 3 

Region 4 

Region 5 

Region 6 

Columns 

3-2 

Gauss, 

e 

g 

- 

1.1430 

0.3780 


0.0730 

0.7630 

- 

0. 3030 

0.7480 

0.6693 


Sphere., Original; 

e 


17.1903 

21.0556 


17.8317 

12.6009 


126.0356 

36.5961 

55.4374 


Sphere, Stable 

wg 

— 

4.3780 

2.2977 


0.9263 

- 0.2114 

- 

30.4585 

11, 7397 

13.4839 


Plate, Original; 

e 


18.6985 

23.9852 


21.1519 

13.9772 


188.8153 

47.0500 

81.0685 


Plate, Stable; 

pWg 

wg 


6.1246 

4.2833 

— 

1.4500 

1.4032 

— 

57.9194 

17.5682 

24.9106 

3-3 

Systematic 

e 

8 

- 

0.3 

0.3 

- 

0.3 

- 0.3 

- 

0.3 

- 0.3 

0.3 


Sphere, Original; 

C 

tiTQ 

— 

0.0882 

0.5409 


0.0094 

- 0.4641 


0.5549 

- 0.4432 

0.4123 

l 

CO 

Sphere, Stable; 

^ws 

ws 


0.1054 

0.8243 


0.2928 

“ 0.6576 


1.9990 

- 0.8188 

0.9895 

CO 

Plate, Original; 

e 

wa 

e 

ws 

- 

0.2552 

0.5636 

— 

0.1789 

- 0.5007 


0.3662 

- 0.5109 

0.4205 


Plate, Stable; 


0.1036 

0.8011 


0.0586 

- 0.6523 


1.9446 

- 0. 8664 

0.9673 

3-4 

Gauss 

e 

e® 

s 

e 

gs 


1.1430 

0.3780 


0.0730 

0.7630 


0.3030 

0.7480 

0.6693 


Systematic 


0.3 

0.3 


0.3 

0.3 


0.3 

0.3 

0.3 


Total 

— 

0.8430 

0.0780 


0.3730 

1.0630 

- 

0.0030 

1.0480 

0.7170 


Sphere, Original; 

gWgS 

wgs 

- 

17.1021 

- 21.5965 

— 

17.8411 

13.0649 


126.5905 

- 37.0394 

55. 7449 


Sphere, Stable; 


4.1861 

2.6868 


1.0686 

0.1490 

- 

30.2402 

11.9819 

13.4409 


Plate, Original; 

^wgs 

wgs 


18.4433 

24.5488 

— 

20.9730 

14.4779 


189.1815 

- 47.5609 

81.2855 


Plate, Stable; 

— 

5.8109 

4.7553 

- 

1.1795 

1.8456 

- 

57.6776 

17.9422 

24.8656 



For 


matrices has upon the magnitudes of the uncertainties in Wg* 

Instance the rms values for e (last column, second and fourth rows) were 

wg 

2 2 

55.4374 W/m and 81.0685 W/m when was computed with the original 

matrices for the sphere and plate, respectively. However, when the 

corresponding stable matrices were used, the values of e were 13.4839 

wg 

2 2 

w/m and 24.9106 W/m , respectively, as seen from the results In the third 
and fifth rows of the last column. 

The next question that arises Is: how do the condition numbers Of the 

stable matrices compare with those of the two original matrices? In TABLE 
4-6, the two condition numbers for each of the four matrices are listed. 

TABLE 4-6. Condition numbers of the original and stabilized matrices for 
the sphere and the plate. 


TYPE OF SATELLITE 

TYPE OF MATRIX 

COMPUTER PROGRAM 

^1 


Sphere 

Original (Unsta- 
ble) 

TARA 7-15-1 

, 131.6 

693.9 

Sphere 

Stabilized 

TARA 7-15-1 

59.9 

223.4 

Sphere 

Partially Stable 

TARA 7-15-2 

71.4 

293.6 

Plate 

Original (Unsta- 
ble) 

TARA 7-15-1 

126.4 

684.7 

Plate 

Stabilized 

TARA 7-15-1 

39.0 

218.2 

Plate 

Partially Stable 

TARA 7-15-2 

58.6 

317.1 


In the first and second rows of TABLE 4-6 are presented the condition 
numbers for the original and stable matrix of the spherical radiometer. 

The fourth and fifth rows show the corresponding condition ntimbers for the 
horizontal flat plate radiometer. These are the four matrices whose results 
have been discussed up to now and which were used In computer program 
TARA 7-15-1. However, the third and sixth rows of TABLE 4-6 present the 
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condition ntindiers of the two matrices obtained by using slightly different 
lower bounds (0.04 for the sphere and 0.02 for the plate) as was discussed 
previously. The results obtained by using these matrices did not meet the 
accuracy requirements and are termed "partially stable" In TABLE 4-6, and 
were of no further use. These matrices were generated by, and used In, 
computer program TARA 7-15-2. 

Data Quality Prediction 

TABLE 4-5 Indicates that the errors In for some of the regions are 
much more significant than those for other regions. For instance, region 5 
shows consistently larger errors than the remaining regions for both 
types of satellites and for all matrices. Hence, even though the results 
obtained with the stable matrix (third and fifth rows of each data group) 
are much more acceptable than those obtained with the original matrices 
(second and fourth rows of each data group)., not all of those results obtained 
with the stable matrices appear equally acceptable. At this point, two 
Important pertinent questions need to be answered. 

(a) What are the accuracy requirements for W^? 

(b) Does the form of the original matrix and the magnitude of its 
elements bear any relation to the errors In W^ obtained for each 
of the regions? If this relationship exists, can It be linked to 
the requirements in (a)? 

These two questions are thoroughly dealt with In the following two sub- 
sections. 

Accuracy requirements .- TABLE 4-7 lists the desired and minimum useful 
tolerances for each of the quantities to be measured for radiation budget 
determinations. It Is seen from this table that the minimum useful 
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TABLE 4-7. Accuracy Requirements 

for Radiation 

Budget Components 

(Recommendations of Investigators conference, 1975) 

Variable 

Accurac 

z 


Desired 

Minimum 

Useful 

Solar Intensity 

± 1 to 2 W/m^ 

± 5 W/m^ 

Solar Spectrum In 

Ozone Bands AX » 50A) 


± 10% 

Components For Global Net: 

Albedo 

± 0.004 

± 0.02 

Longi/ave Exltance 

± 1 W/m^ 

± 5 M/m^ 

Components For Regional Net: 

1 

Albedo 

± 0.02 

± 0.08 

Longwave Exltance 

± 3 W/m^ 

± 15 W/m^ 

Medium Resolution Scanning: 

Albedo 

± 0.04 


Longwave Exltance 


± 6 W/m^ 



Long Term With 
Monthly Resolution 
Desired - Seasonal Is 
Minimum Useful Period 


Monthly For 10° Of Great 
Circle Latitude And 
Longitude 


Monthly Averaged Determined 
From Scanning Data - 10^ 
to lOS km2 Spatial Resolution 



accuracy for monthly averages of the longwave exitance (W^) of small regions 
2 

Is ± 15 w/m . In the cases presented In TABLE 4-5, however, the V£dues of 
are Instantaneous for all practical purposes since the time Interval 
during which the measurements were taken is of the order of minutes. Hence, 
the tolerances for these results should be less stringent than those for 
monthly averages stated above.. 

It should be pointed out that the accuracy requirements listed in 
TABLE 4-7 are the results of recommendations made at the Chicago Investi- 
gators conference of 1975. On the basis of the above tolerance of ± 15 
2 

w/m , one can see that from all the values obtained with the two 
stabilized matrices the only ones which are not acceptable are 

(a) Those of region 5 resulting from gausslan, or combinations of 
gausslan/systematlc power errors for the spherical satellite. 

(b) Those of regions 5 and 6 resulting also from gausslan or 
gausslan/systematlc power errors for the plate. 

The possibility of a connection existing between the structure of 
the original configuration factor matrix and the stability of its inverse 
matrix was thoroughly Investigated and excellent results were obtained as 
discussed In the following paragraphs . 

Matrix Parameters .- The original configuration factor matrices shotm in 
figures 4-3 and 4-4 for the sphere and plate, respectively, are again intro- 
duced here in figures 4-9 and 4-10. However, these latter figures display 
additional Information which will be needed In this discussion. Each of 
the first six quantities In the last column consist of the sums of the 
elements in their corresponding rows, and should be approximately equal to 
the shape factor of the total FOV of the observation corresponding to that 
row. The seventh quantity In the same column Is the sum of the first six 
quantities. When this sum Is divided by six, it yields the average shape 
factor for the total FOV, which is shown in the last row of the last 
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OOfSPH(K,p OOSSPH(K,2) OOFSPH(K,3) 

0.484847428 0.384856899 0.127704158 

0.276829292 0.276829292 0.276829292 

0.108360067 0.127704158 0.384856899 

0.017068977 0.035500329 0.355753027 

0.000000000 0.000000000 0.247608706 

I 

I 

0.000000000 0.000000000 0.137012641 

COLUMN SUMS 0.887105765 0.824890678 1.529764721 

PERCENTAGES,'X 81.416314223 75.706371568 140.39792|l2JL2 


Figure 4-9. Original matrix for the spherical radiometer 


COFSPH(K,4) 

C0FSPH(K^5) 

G0FSPH(K 6) 

SPHROW(K) 

0.108360067 

0.000000000 

0.000000000 

1.105768551 

0.276829292 

0.000000000 

0.000000000 

1.107317169 

0.484847428 

0.000000000 

0.000000000 

1.105768551 

0.595573384 

0.017068977 

0.035500329 

1.056465024 

0.596414856 

0.068729515 

0.143711947 

1.056465024 

0.415871635 

0.137012641 

0.415871635 

1.105768551 

2.477896661 

0.. 222811133 

0.595083912 

6.537552870 

227.415062849 

20.449039970 

54.615290183 

1.089592145 
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COFHOR(Kil) 

CDFH0R(K,2) 

C0FH0R(K 


C0FH0R(K,4) 

C0FH0R(K,5) 

C0FH0R(Kj6) 

BORROW (K) 


0.378000655 

0.279486438 

0.078892^73 

0.063693275 

0.000000000 

0.000000000 

0.800073041 


0.200205230 

0.20020523 

0.200205230 

0.200205230 

0.000000000 

0.000000000 

0.800820918 


0.063693275 

0.078892673 

0.279486438 

0.378000655 

0.000000000 

0.000000000 

0.800073041 


0.008805215 

0.018591209 

0.251441001 

0.471341713 

0.008805215 

0.018591209 

0.777575561 


0.000000000 

0.000000000 

0.172299495 

0.473440034 

0.040187252 

0.091648780 

0.777575561 


0.000000000 

0.000000000 

0.087326867 

0.312709653 

0.087326867 

0.312709653 

0.800073041 

COLUMN SUMS 

0.650704375 

0.577175549 

1.069651704 

1.899390561 

0.136319334 

0.422949643 

4.756191165 

PERCENTAGES. Y 

82.087244037 

72.8114824591 

34.9380207892 

39.610708856 

17.196869788 

53.355674072 

0.792698527 


Figure 4-10. Original matrix for the horizontal plate radiometer. 

■ i 



colunm. This average value Is denoted by x for the sphere, and by y for the 
plate. Each of the quantities In the row entitled "column sums" is the 
sum of the six quantities above it. The first six elements of the row termed 
"percentages, x,y" are obtained by dividing each colxmn sum by the average 
shppe factor for the total FOV and multiplying the result by 100. 


NOTE: Since the systematic errors do not produce unacceptable results, 

the following discussion applies only to those cases in which the 
errors are totally or partially gausslan. 


Prediction scheme .- A qualitative prediction classification based solely 
on the structure of the original matrices was found. It yielded excellent 
results when compared with the xms's of several errors computed for 
each of the regions. This rms is based on ten sets of error data; each set 
has six e values required to perturb the six power measurements of each 
observation group. 


In order to predict the quality of the data to be retrieved, three 
distinct data classes were arbitrarily selected. The quality of the 
values to be retrieved were predicted to be either ACCEPTABLE; POOR, or 
REJECTABLE, according to the criteria detailed below. 

Let SS^ denote the sum of the elements in the Jth column of the 
configuration factor matrix for the sphere, and let SP^ denote the equivalent 
sum for the plate. It is recalled that x was used to denote the average of 
the six shape factors for the total FOV of the sphere. Similarly, y was used 
to denote the equivalent quantity for the plate. Then, using these 
definitions, the first criterion can be stated as follows: 


NOTE: In the following, REJECT means reject the value determined for the 

jth region from observations by the sphere, in the case of SSj , or by 
the plate, in the case of SP^ . Similar meanings should be attached 
to ACCEPT and POOR. 
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IF SS. < 0.2X 
J 


REJECT 


(4-53) 


IF SPj < 0.2Y 


REJECT 


Otherwise, SS^ and SP^ 


are subjected to the following tests. 


(4-53’) 


IF SS, > 1.25X 
J 


IF SP, > 1.25Y 
3 


ACCEPT 


ACCEPT 


(4-54) 

(4-54’) 


If this test is not passed, the diagonal element F^^ of the j-column of both 
matrices (F^^ for the sphere and F^^ for the plate) are. subjected to the 
following test. 


IF F.. <0.25 SS. 
Jj J 


REJECT 


(4-55) 


IF F“. < 0.25 SP. 
33 - 3 


REJECT 


(4-55’ ) 


Otherwise, the two sets of quantities are subjected to a final test. 


IF F.. > 0.6 SS. 
jj 3 


ACCEPT 


(4-56) 


IF F?. > 0.6 SP. 
33 3 


ACCEPT 


(4-56’) 


For those cases when Fj. and F^ are between the above two limits, the 


jj jj 


classification is. 


IF 0.25 SS. < F. . i SS. 

3 33 3 


POOR 


(4-57) 
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POOR 


(4-57») 


IF 0.25 SP. < F?. < 0.6 SP, 

3 22 2 

The predictions made on the basis of the above criteria were sufficiently 
satisfactory to indicate that there is a definite relationship between the 
quality of the data to be retrieved and the configuration factor matrix of 
the measurements. This is of great value since a simple analysis of the 
matrix tells in advance which are the regions whose results should be 
considered for further processing, as will be discussed in the following 
subsection. 

TABLE 4-8 compares the data quality predictions with the rms's of the 
errors in for each of the regions. As seen from this table, in every 
instance the predictions agree with the computed errors in W^. From 
TABLE 4-8 the following can also be seen, 

(a) In all cases the uncertainties in for the plate are larger 

than those for the sphere 

(b) The only Instance in which a REJECTION was predicted was for the 

plate in region No. 5. The errors computed for the plate in 

2 

this case were, 33.9284 w/m for the gaussian perturbation and 
2 

33.8025 w/m for the gauss Ian/ systematic error combination. 

It is apparent then that before retrieving the values of from a given 
set of power measurements, the above scheme will tell which values will 
be good enough for further processing, such as averaging of data discussed 
in the next subsection. 


Weighted Averages 

The n observations required to solve the n simultaneous equations for 
can be taken within a time interval of any arbitrary length. However, 
the physical situations represented by the results would be different for 
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TABLE 4-8. Comparisons of data quality predictions with rras’s of computed errors. 


Errors (e .) In W (W/m^) 

w C 

1 Region 2 Region 3 Region 4 Region 5 Region 6 


Type of Type of 
Satellite Power Error 
(W/m2) 


Type 

of 

Results 


Region 


Sphere 


Gaussian, 


Predicted Poor 


Sphere 

Gaussian, 

e 

g 

Computed 

15.2497 

Sphere 

Gaussian Plus 
(0.9 W/n»2) 
Systematic, 

e 

ga 

Computed 

15.1991 

Plate 

Gaussian, 

e 

Predicted 

Poor 



g 



Plate 

Gaussian, 

c 

g 

Computed 

15.5824 

Plate 

Gaussian Plus 
(0.9 W/m2) 
Systematic, 

e 

ga 

Computed 

15.4480 


Poor 

Accept 

Accept 

Poor 

Accept 

19.9081 

6.8638 

2.9748 

17.4476 

5.7264 

20.0824 

6.8882 

/■ 

3.1020 

17.4142 

5.8314 

Poor 

Accept 

Accept 

Reject 

Accept 

21.6890 

9.0607 


33.9284 

9.0240 

21.9173 

9.0688 

4.2407 

33.8025 

9.2882 



different time Interval lengths. At any rate, regardless of the magnitude 
of the time Intervals selected for accomplishing each set of observations, 
the resulting data can always be averaged out over much longer time periods. 
For purposes of the present discussion, a set of results will be considered 
Instantaneous If the corresponding set of n observations Is taken during 
a single pass or orbit of the satellite. A total of X sets of data will 
be considered for averaging In the following discussion. 

The following three types of averages were considered. 


W . 
ei 


Weij 


T 


(4-58) 



ill ^eij 
T 

jSl Aij 


(4-59) 


ex 


-jil ^eij 
T 

jii 


(4-60) 


Where, 



the jth determination of the value of for the 1th 
region 



“ total area seen of the ith region during the jth 
observation set. 


F 

T 

W 


ij 


el 


a sum of all the configuration factors of the 1th 
region which entered in the jth set of observations. 

a total number of results, or observations sets, used 
in the averaging. 

a plain average of the value of W^ for the ith region. 
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W^.“ Weighted average of W for the ith region, using the area 
as weight. 

F 

Wei* Weighted average of W for the ith region, using the 
^ configuration factor as weight. 

The average given by (4-58) , assigns equal weights to all the W^ 
results entering into the averaging process, regardless of the sizes of 
the segments of the ith region that were within the FOV's of the radiometer, 
and regardless of the positions that these segments occupied within the 
radiometer's FOV. Hence, this type of averaging was considered inadequate. 

The type of average defined by (4-59) takes into account the size of the 
region segments that entered into the observations; however, two segments 
of equal area but making different contributions due to their different 
positions within the FOV are given equal weights, which might be undesirable 
in some cases. The third type, defined by (4-60) was considered adequate since 

it does_not-Suffer_from either-of— the-shortcomings mentloned-above— for 

(4-58) and (4-59). 

Thus, weekly or monthly weighted averages (for any time of the day) 
can be easily obtained from Instantaneous results (for the time of day 
selected) by using the expression (4-60) . Although the SWR component 
is not discussed in detail in this report, a proposed method that is 
currently under investigation for tackling the problem of reflected SWR 
will now be discussed. 

Proposed Procedure for Computing W^ 

When the W^ regions of the E-A system previously defined were consi- 
dered, it was assumed that the variations of W^ within the region were 
negligible or barely detectable by the satellite radiation sensing 
system. Now, one can assume another common characteristic to all 
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area elements within a given region; this Is the ALBEDO ACC^) for a specific 
solar zenith angle For reasons which will be apparent later, the 

zenith angle chosenfbr this common characteristic Is » 0°. Before pro- 
ceeding further. It Is advisable to recall the following relationships that 
were Introduced In the section entitled "BASIC IDEAS". 

W ( 5 ) = difi N (9,i(),5) sin 0 cos 0 d0 (2-19) 

r 0 0 ^ 

A(C) = r(C) = di|i p(d,ip,z) sin 0 cos 0 d0 (2-20) 

0 0 


A(5) = r(0 = W^(0/H(5) 


( 2 - 21 ) 


p(0,i|;,O = N (0,<|»,?)/H(?) 


( 2 - 22 ) 


Then, one can write for those characteristics that have been assumed 


common to all area elements within a given region the following. 


W (0)= dij> N (0,iJ/,O)sin 0 cos 0 d0= constant (4-61) 

^ 0 0 ^ 

A(0)- r(0)= di|) p(0,ij;,O)sin 0 cos 0 d0= constant (4-62) 

o 0 

However, since at any Instant of time the solar zenith angle exhibits 
spatial variations within any given region; then, r(c), A(^), and W(^) are 
not necessarily constant within the given region. Then, the problem is, how 
is one to determine the W^(q) of a region If the value of W^( 0 ) for that 
region Is known? 

The following Is a simple procedure which uses data available from 
previous satellite observations to provide a quick solution to the above 
problem. 
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Figure 5 on page 11 of reference 4 shows several curves which depict 
relationships between r($) and r(6) for different types of surface regions 
of the E-A system. One or several of these curves are selected to repre- 
sent the reflecting characteristics of the hypothetical regions one has 
assumed. Consider now the 1th area element of the kth region which Is 
within the FOV of the jth observation of a satellite radiometer. This area 
element is denoted by and the radiant power it contributes to the 

jth measurement Is The reading from the curve in figure 5 that 

corresponds to is denoted by 




(4-63) 


The albedo A and directional reflectance r for 4A^^j^ is given by 




(4-64) 


It should be noted that i^j^(u) does not have the subscripts 1 and j 
since the value rj^(o) " changes only if the region changes. 

From (2-21) shown above, one writes 




( 4 - 65 ) 


Then, substituting (4-64) into (4-65) one obtains 




(4-66) 


The expression for the power increment "which the area element 

emits in the direction of the satellite is. 


^I’ljk 


\ijk(®»’l'»^ 2 


A A. cos 6 . 
X]k i3k 


(4-67) 
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where. 




Reflection radiance of AA... In the direction 
(6, \l>) due to solar radlaHon Incident from a 
direction given by the solar zenith angle C. 


A » ira 
8 


Characteristic area of the radiometer, where a 
Is the radius of the radiometer. 


ijk 


Distance from radiometer. 


ijk 


Zenith angle of the radiometer as seen from 


From (2-22) shown above, one writes 




(4-68) 


Substituting (4-68) into (4-67) and assuming that A > Im , one obtains 

s 




A A, .. cos 9 , 
Ilk nk 

4)k 


(4-69) 


Now, from the corresponding curve in one of the figures B-3, B-4, or B-5 
of reference 3, one can read directly (or by using linear extrapolation) the 
value of the ratio r(c)/Tr p(0, C) which will be denoted by R’ , that is. 




(4-70) 


or. 






(4-71) 


Substituting (4-71) Into (4-69) , one obtains 


= H,,, (O 




Ijk “ijk^'’' TT R^^j^(9,ij),0 ^2 


A A. cos 6 . 

1.1 k nk 


ijk 


(4-72) 
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But from (2-21) above, one writes 




Substituting (4-73) into (4-72) one obtains 


(4-73) 


A A. 


AP 


ijk 


ijk 


). C 


cos 6 


ilk 






(4-74) 


By the definition of the shape, or configuration factor given. ^ by (2-2), 
2 

with A “ Im , one writes 
s 


AP. = F, W . „ (c) 
ijk ijk rijk 


Comparison of (4-74) and (4-75) indicates that 


(4-75) 


AA. 

F..,. - ( 


COS 9 . . , 

iiL 


) 13 J 

Rl„ (0,'P,?) dj,, 


(4-76) 


ijk V 

''ijk"'”'”' “ijk 

In the following discussion, a new type of shape factor will be Introduced 


in order to develop a simple expression for computing the solutions for the 
reflected SWR component of the energy budget. This new type of shape factor 
will be termed PSEUDOSHAPE FACTOR. 

One can represent H (O as 


= \ W «ljk 


(4-77) 


Note that since the value of Rjj.(o) changes only with changes of region, 
only Che subscript k of the region need be used, just as in the case of 
rj^(o) and Aj^(o). 


4-49 



Substituting (4-77) into (4-65) one obtains 


“r<°> ‘ijk 

For the case ^>0, one obtains from this expression, 


(4-78) 




(4-79) 


Again here, only the subscript k is necessary since these three quantities 
change value only with region changes. 

From (4-63) one writes 


'ljk"> + %4k<'> ■^k<« 

Substituting (4-80) into (4-78) one obtains 


(4-80) 


"rijk<^> 'k<0' «k«» 

And substituting (4-79) into (4-81) , 


(4-81) 


"rijk<^> ^ “rk<°) \jk<^) «lik 


(4-82) 


Substitution of (4-82) into (4-74) yields 


AA R. (5) cos 9,., cos 

IP, . ( -tis. ): uiL , (0) 


IJk 


7T 


Rijk <* Ijk 


rk' 


From this expression one defines the PSEUDOSHAPE FACTOR 


(4-83) 


AA... R. ., (?) cos 6... cos 

F!. -( 


ijk 


TT 


^ijk ‘^ijk 


(4-84) 
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By conparlng the expressions (4-76) for (4-84) for o®® sees 

that these two qualities are related by the following expression. 


F, = F.., R.., ( cos c- 
xjk xjk xjk ^xjk 


(4-85) 


Let the following quantity be defined. 


^ijk 

Then, (4-85) can be' rewritten as 


(4-86) 


F! = F. Rl (c) 
xjk xjk xjk ^ 


(4-87) 


For future use, the weighted average value of (^) over all those 

area elements 4A. , of the hth- region which appeared at least once within 
ijR 

one of the POV*s of the radiometer will now be calculated. Thus, one 
writes 


I J 


E E F' = E F... R!I, (O = R," E F. 

i=l j=l ilk i,j ijR ijR ^ i,j 


(4-88) 


where. 




total number of area elements of the kth, region 
within the FOV of the jth observation. 

the number of observations in which the kth region 
appears. 

weighted average of R!!> ( 5 ) for the kth region. 

IJ X 


Thus, 



Returning to equation (4-83) , one obtains the power contributed by the 
kth region to the jth observation, by adding up over all the area elements 
In the kth region. Hence, 

I 

P., = E 4P. 

jk i=i ijk (4-91) 


Then, one obtains the total power detected In the jth observation by adding 
up all the powers P^j^ contributed by the K regions, that Is 


K 

P. » E P., 
J k=l jk 


(4-92) 


These Pj'^ elements of the power column matrix {P}. In exactly the 

same manner, one obtains for the configuration factor and the pseudo- 
configuration factor that the kth region contributed to the jth 
obseivatlon 



I 

1=1 Ijk 


(4-93) 



I 

E f! 
1=1 13 k 


(4-94) 


These constitute the elements of the pseudoconfiguration factor 

matrix F'. Thus, one can write in matrix form. 


F'{W^(0)} = {P} (4-95) 

where the elements of the column matrix {W^(o)} are the hypothetical values 
of W^(o) which were assigned originally to the different regions. Then, 
by operating on {P} with F'- the inverse of F' , one obtains back (W^(o)}, 
that Is 
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(4-96) 


F'"^{P} = {W^(0)> 

If this matrix operation retrieves the correct values of W^(o) that 
were originally assigned to each of the regions, then one proceeds, as In the 
one of W , to perturb {P} with gausslan and systematic errors to test the 
stability of F 

The next step Is to obtain the values of W^(c) for each of the regions 
from the corresponding values W^(o) found through (4-95). This can be done 
as follows. 

Assisnlng that the hth. region appeared In the FOV's of J observations, 
one can find the total power that this region contributed to the J 
observations by adding over j the P^j^ powers given by (4-91) , that is. 


J 

Fr = 2 P„ 

k j=l 3k 


(4-97) 


Similarly, the contribution of all the configuration and pseudoconfiguration 
factors by the kth region to all the J observations are, from (4-93) and 
(4-94) , 


^k = fii ^jk 


(4-98) 


F» = E F! 
k j=i jr 


(4-99) 


And from (4-91) and (4-97) is obtained 


P, = 2 AP.., 

k ijk 

By comparing this expression with (4-75), one can write. 


(4-100) 


P 


k 


^ - 1 , W 

j,i ijk rxjk 
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(?) 


(4-101) 



If one defines to be the weighted average of W^(5) over the kth 
region for all the J observations, one writes 




rk 


.2. Fijk 


(4-102) 


Then, (4-101) becomes 


P, » W , Z F. 
k ijk 


(4-103) 


If now one compares expression (4-100) with equation (4-83) and uses (4-84) , 
one can write, 


P, = W , (0) Z F!„ 
^ j,l ijk 


Equating (4-103) and (4-104) , one obtains 


(4-104) 


W , = 
rk 


Z FI 

2 F,- -1, 

i.j 


“rkW> 


(4-105) 


But the ratio on the right Is the same as that In (4-90). Hence, one can 
say that 


“rk = “rk<°) K' W-10« 

Where Is the weighted average of r( ^ cost. /r(o) for the kth region for 
all observations. 

Therefore, the essence of the procedure when applied to an actual set 
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of power measurements is as follows. Compute the elements of the configu- 
ration and pseudoconfiguration factor matrices. Obtain the inverse matrix 
of the pseudoconfiguration factor matrix and operate with it on the 
column matrix of the power measurements in order to obtain the column matrix 
{W^(o)}. From these elements, the corresponding elements are obtained 
by application of (4-106) . 
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BEST PIT/INVERSION TECHNIQUE 


As was mentioned in the last paragraph of the section entitled "FUNDA-_ 
MENTIS OF THE TECHNIQUES the main feature of the Best Fit/Inversion 
Technique is that the number of observations m is larger than the number 
of unknowns n. An extension of the method of least squares is used to 

determine the - approxlmajtingj surface in the nrhl dimensional space.' The 
procedure followed to find the values of is Illustrated by a simpli^; 
fled three dimensional case as follows. 

Let tt. be assumed that a satellite has made eighteen observations over 
two regions. The POV of each observation is either totally filled by one of 
the reglonsy or totally filled by segments of both of the regions. Further- 
more. let it be assumed that the W values of the two regions, identified as 

B 

regions 1 and 2, are in the ranges » 240.0 ± 5 w/m and W ^2 ® 280.0 ±5 
2 

w/m , as shown^ in TABLE' 5-1. This table shows the configuration factors 
^1 ^2 regions 1 and 2 contribute to each of the measurements. The 

hypothetical values of and W ^2 the time each observation is taken are 
also presented in this table, as well as the two partial powers for each of 
the measurements. These partial powers, for the ith observation, are given by 


P. = P. -+ P. W . - F, , A +W . .F. .A 
X x,l x,2 ei,l i,l s ex, 2 x,2 s 


(5-1) 


Where 


A; « Im ■ characteristic area of the radiometer, 
s 


Therefore, the eighteen observations can be written as 

el,l Ipl el , 2 lp2 1 


”e2,l ^2,1 ^2,2 ^2,2 “ ^2 


(5-2) 


W F + W F = p 
e3,l 3,1 ^ e3,2 *’3,2 ^3 


^18, l'*' \l8,2 ^18,2 “ ^18 
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TABLE 5-1. Hypothetical data about 

satellite 

observations of 

two regions. 





Measurement No. 










Parameter 

1 

2 

3 

4 

5 

6 

7 

8 

9 


O.A 

0.5 

0.6 

0.3 

0.7 

0.2 

0.8 

0.1 

0.9 

„^2 

0.6 

0.5 

0.4 

0.7 

0.3 

0.8 

0.2 

0,9 

0.1 

F « ?! + Fj 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

( 
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Wei ^ 

240 

241 

242 

243 

239 

238 

237 

244 

^ 245 

W^2 

280 

279 

278 

277 

281 

282 

283 

276 

275 


(w) 

96.0 

120.5 

145.2 

72.9 

167.3 

47.6 

189.6 

24.4 

220.5 

?2 (w) 

168.0 

139.5 

115.2 

193.9 

84.3 

225.6 

.56.6 

248.4 

27.5 

P “ Pi + P2 

264.0 

260.0 

260.4 

266.8 

251.6 

273.2 

246.2 

272.8 

248.0 



TABLE 5-1. 

Hypothetical data about 

satellite 

observations 

of two regions. 

(Continuation) 



Measurement 

Parameter 

No. 

10 

11 

12 

13 

14 

15 

16 

17 

18 



0.0 

1.0 

0.45 

0.55 

0.65 

0.25 

0.35 

0.75 

0.15 

X. 

F 

^2 


1.0 

0.0 

0.55 

0.45 

0.35 

0.75 

0.65 

0.25 

0.85 

F « Fi + 


1.0 

1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


Ut 

1 

Oj 

^el 

(w/m^) 


240 

236 

235 

240 

241 

242 

243 

239 

238 

”e2 

(w/m^) 


280 

284 

285 

280 

281 

282 

283 

279 

278 


Pj (w) 

0.0 

236.0 

105.75 

132.0 

156.65 

60.5 

85.05 

179.25 

35.7 

P 2 (w) 

280.0 

0.0 

156.75 

126.0 

98.65 

211.5 

183.95 

69.75 

236.3 

P - P^ + 

280.0 

236.0 

262.50 

258.0 

255.30 

272.0 

269.0 

249.0 

272.0 



The 1th equation of this set represents the 1th plane In a three-dimensional 
space whose rectangular coordinates are ^1* form of the 

equations Indicates that all planes cross the origin of the coordinate system. 

One asstimes further that the data obtained from the observations (l.e., the 
values of ^1^!2 ^1^ when plotted In this three dimensional coordinate 

system become a scatter diagram which can be represented by an APPROXIMATING 
PLANE (see reference 9) that best fits bhe data. The equations of this plane 
will be of the form 


P 





W 


e2 2 


(5-3) 


and the partial regression coefficients and will be determined by an 
extension of the method of least squares as follows. 

The general form for the equations (5-2) ls» 


P, = W , - F, , + W . „ F. - 
1 el,l 1,1 ex, 2 x,2 (5—4) 

One multiplies this equation through by F. . to form 18 equations (1^1,... 18) 
These 18 equations are added up to generate the first equation of the set of 2 
equations needed to find the 2 unknowns. The second equation required Is 
obtained by multiplying (5-4) by F^_2 to form, again, 18 equations which are then 
added together as before. Thus 


18 



18 

2 

18 

E 

P,F. , 

= w 

E 

FT W ^ 

5^ F , F „ 

1=1 

1 1,1 

el 

1=1 

i,l e2 

i=l 1,1 "^1,2 

18 

P.F. ^ 

= W 

18 

E 

F. , F. ^ 

18 2 

+ W „ E F „ 

I 

1=1 

1 1,2 

el 

1=1 

x,l 1,2 

e2 x,2 


(5-5) 

(5-6) 


In order to simplify the notation, the following quantities are defined. 
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_ 18 

PF » Z P. 1 
1 i=i 1 i,l 


18 

PF„ = Z P. F. 


i=l 


i i.2 


2 18 2 

f: = Z FT . 

1 i=i x,l 

f2 « Z F? „ 

2 i=i x,2 

18 


F F» = ^ F F » 

r2 i=i i,i ^i,2 

Using (5-7) , equations (5-5) and (5-6) can be written aiS 


(5-7) 


PF, = W 1 F- + W - F F- 
1 el 1 e2 1 2 


PF =W FF + W F 
2 el 12 e2 2 


(5-8) 


As was mentioned before, (5-8) is a set of two equations in the two 

unknowns and which is required to solve the problem. 

A short program was written to solve the set of equations (5-8) by using a 

matrix inversion subroutine. This program also computes the average values of 

and W ^2 order to compare these results with the coefficient values found 

by the least squares method. As in previous programs, the original matrix and 

its inverse were multiplied in order to verify that the result was the identity 

matrix. The values of W , and W « obtained from the eighteen observations are 

el e2 

compared in TABLE 5-2 with the average values of the hypothetical distributions 

of U . and W . given in the problem, 
el ez 


TABLE 5-2. Comparison of values computed from error- free data with the 
averages of the given values. 


TYPE OF VALUES 

W^^(wZm^) 


By least square method from 

239.83 

279.99 

18 observations 



Averages of the given 

240.17 

280.17 

values' 
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Tvo important modifications were introduced in the next application. 

(a) Six regions, rather than two, were considered and a total of thirty -six 
observations were made, (b) systematic and gausslan uncertainties were 
Incorporated into the power measurements. It should be pointed out...that.. since 
the space now considered is 7th dimensional, one can no longer speak of the 
equation of an approximating plane, but rather of the equation of an 
approximating surface. 

2 

A gausslan random error distribution, with a sigma value of 0.5 w/m. 

2 2 

was used, as well as three values of systematic errors, 0.3 w/m , 0.6 w/m 
2 

and 0.9 w/m . 

TABLE 5-3 lists five groups of results based on the following types of 
uncertainties in the observations. 

1. Ho error in the observations. 

2. Gausslan errors only included. 

2 

3. Systematic (0.3 w/m ) combined with gausslan errors. 

2 

A. Systematic (0.6 w/m ) combined with gausslan errors. 

2 

5. Systematic (0.9 w/m ) combined with gausslan errors. 

The five sets of data in TABLE 5-3 clearly show", that the results are 
not highly sensitive to the gausslan and systematic uncertainties assumed for 
the power measurement. 

It is recalled that the essential difference between this technique and 
the Instantaneous/Inverse Technique is that in the latter the number of 
observations m equals the number of unknowns n; while in the former, m is 
greater than n. However, the basic equations to be used are the same in both 
cases, namely, equations (3-22) and (3^23) introduced in the section entitled 
"Mathematical Bases." The configuration factors in these two expressions 
include the LDF as can be seen from equations (3-18) and (3-19). Furthermore, 
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TABLE 5**3. Comparison of values computed from data having uncertalntltles with the averages of the 
given values. 


VALUES OF Wg (w/m^) 


TYPE 

OF DATA 

Region 1 

Region 2 

Region 3 

Region A 

Region 5 

Region 6 

AVERAGE of Given Values 

179.36 

200. A2 

220. AA 

239.83 

259.56 

280. A A 

FROM 

Error-Free Data 

178.09 

199.63 

221.22 

2A2.15 

256.66 

280.69 

FROM 

Observations with 2 

Gaussian Errors (o«0.5 w/m ) 

178.38 

199.65 

221.30 

2A2.06 

256. A5 

280. 6A 

FROM 

2 

Observations with (0.3 w/m ) 
Systematic & Gaussian Errors 

178.68 

199.95 

221.60 

2A2.36 

256.75 

280. 9A 

FROM 

2 

Observations with (0.6 w/m ) 
Systematic & Gaussian Errors 

178.98 

200.25 

221.90 

2A2.66 

257.05 

281. 2A 

FROM 

2 

Observations with (0.9 w/m ) 
Systematic & Gaussian Errors 

179.28 

200.55 

222.20 

2A2.96 

257.35 

281. 5A 



the corresponding expressions for the case of reflected SWR are developed In 
a manner similar to that employed In the Instantaneous/Inverse Technique. 

The results obtained so far by this technique are satisfactory and further 
work In Its application should be pursued. 
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SPHERICAL EARTH-ATMOSPHERE SYSTEM 


Once it had been shown that the techniques for computing yield 
acceptable results. It was decided to return to the original problem In 
which a spherical earth-atmosphere system was being considered. 

As was mentioned previously, the power Integrations are to be accomp- 
lished by adding up the power increments reaching the satellite simultane- 
ously from the different area elements whlthln the FOV. Hence, the first 
step was to generate a computer program which would divide the surface area 
of the earth-atmosphere system Into a finite number of area elements of 
equal area. The data output of this program gives the longitudes and lati- 
tudes of the centroid and four boundaries of each of the 2060 area elements. 

2 

An area element of 250,000 km was found to be adequate* This size corre;- 

sponds to an area of 500 by 500 km, or approximately 5° by 5° of great circle 

arc. This yielded a total of 2060 area elements, 2058 of which have equal 
11 2 

areas (2*5x10 m each) . The other two elements are those centered at each 

11 2 

of the poles. These have an area of 2.339x10 m each. This yields a total 

14 2 

surface area for the spherical E-A system of 5.149678x10 m . The radius 

of this system is considered to be 6*40155xl0^m, which gives a surface area 
14 2 

of 5.149679x10 m . Since the radius of a sphere having an area equal to 
the area of the earth has a radius of 6.37123xl0^m, it follows that an atmo- 
spheric spherical shell of 30.32 km thicness is being included* This is 
approximately the thickness commonly used (30km) since it was first suggested 
by Dr. Frederick B. House (reference 10). More than 95% of the total atmo- 
sphere is contained within this shell. 

The output data of this program is used in the main computer program 
to calculate the values of a spherical E-A system* 
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The total shape factors for the sphere and plate were calculated also in 
this main program by adding up the shape factors of each of the area 
elements within the FOV of the radiometers. The results are compared 
below with those obtained analytically. 

Configuration factor of the total FOV 
Radiometer Numerically Analytically 

Sphere 1.089592145 1.0838471 

Plate 0.792698527 0.7901660 

These results Indicate that It Is unnecessary to use area elements of 
smaller size to perform the numerical Integrations. 

The results of the Investigations presently being performed using a 
spherical £-*A system will be discussed in a subsequent report. 
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CONCLUSIONS 


On the basis of the investigation results reported in this document, 
the following has been concluded: 

1. The Instantaneous /Inversion Technique (including the data 
quality prediction and matrix stabilization schemes) yields 
excellent results when applied to radiometer data containing 
gausslan and systematic errors. The values obtained are 
acceptable according to the pertinent accuracy requirements 
displayed in TABLE 4-7. 

2. The problem of determining the solar radiation reflected by the 
E-A system has been cast in a form which requires minimum com- 
puter time for calculating and the albedo. It is concluded 
that this simple formulation represents the optimum method to 
solve this problem and will soon be implemented. 

3. The results obtained from preliminary applications of the Best 
Fit/Inversion Technique clearly indicate that the errors in 
calculated by the use of this technique are not highly sensitive 
to gausslan and/or systematic power uncertaintltles . In all- 
cases investigated the values retrieved met the pertinent 
accuracy requirements in TABLE 4-7. 

4. It has been concluded that division of the surface of the E-A 

5 2 

system into 2060 area elements of about 2.5x10 . Vm each yields 
accurate results when used in numercial integrations of 
configuration, factors. 
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